Q3-3--J 


[< 


.  AFCRL  63-429 

Hf 
lv 

oo 
o 

rfi 


kj 

BET* 

QQ 


O 

L,U  ta*«=» 

C E? 

SQ 

J?C  GO 

P  C 


THE  KINETIC -ENERGY  BUDGET  OF  THE  STRATOSPHERE 
OVER  NORTH  AMERICA  DURING  THE  WARMING  OF  1957 


by 

M.  A.  LATEEF 

Department  of  Meteorology 
Florida  State  University 


Hi 

»l 


i 


i 


Contract  No.  AF  19(628)  -394 
Project  No.  8604 
Task  No.  860405 


SCIENTIFIC  REPORT  NO.  1 
15  April  1963 


AST ! A' 

I  8AAY  6  1963  1* 

Jy  i±i  O-Ain  U  V 

4  TIStA 


PREPARED  FOR 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
OFFICE  OF  AEROSPACE  RESEARCH 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS 


AFCRL  63-429 


4 


THE  KINETIC- ENERGY  BUDGET  OP  THE  STRATOSPHERE 
OVER  NORTH  AMERICA  DURING  THE  WARMING  OF  1957 


by 

M.  A.  LATEEF 

Department  of  Meteorology 
Florida  State  University 


Contract  No.  AF  19(628)-394 
Project  No.  8604 
Task  No.  860 405 


SCIENTIFIC  REPORT  NO.  1 
15  April  1963 


Prepared  for 

Air  Force  Cambridge  Research  Laboratories 
Office  of  Aerospace  Research 
United  States  Air  Force 
Bedford,  Massachusetts 


Requests  for  additional  copies  by  Agencies  of  the  Department  of 
Defense,  their  contractors,  and  other  Government  agencies  should 
be  directed  to  the: 

ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 

Department  of  Defense  contractors  must  be  established  for  ASTIA 

services  or  have  their  "need-to-knovM  certified  by  the  cognizant 

military  agency  of  their  project  or  contract. 

All  other  persons  and  organizations  should  apply  to  the: 

U.  S.  DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TECHNICAL  SERVICES 
WASHINGTON  25,  D.  C. 


ACKNOWLEDGMENTS 


The  author  is  grateful  to  Professor  Richard  A.  Craig  for  his 
guidance  throughout  the  course  of  this  investigation.  Thanks  are 
extended  to  Messrs.  W.  J.  Koss,  R.  S.  Hawkins  and  P.  R.  Sfcicksel  for 
aid  and  cooperation  during  the  programming  of  the  IBM  709  computer. 

A  portion  of  the  programming  costs  of  this  investigation  was  home 
by  the  National  Science  Foundation. 

Appreciation  is  also  due  to  Professor  R,  J.  Reed  of  the  Uni¬ 
versity  of  Washington  for  communicating  to  the  author  results  of 
some  of  his  computations. 


ii 


TABLE  OF  CONTENTS 


Page 

ACKNOWLEDGMENTS .  ii 

TABLE  OF  CONTENTS .  iii 

LIST  OF  TABLES  . v 

LIST  OF  ILLUSTRATIONS . .  *  .  .  vii 

ABSTRACT . . .  1 

Section 

1.  INTRODUCTION .  3 

2.  THE  MECHANICAL  ENERGY  EQUATION . .  ..  8 

Derivation 

Physical  meaning  of  the  terms 
Previous  applications 

3.  PREVIOUS  INVESTIGATIONS . .  .  .  1? 

4.  COMPUTATIONAL  PROCEDURES .  22 

Area  and  line  integrals 
Vertical  integration 
Input  data 

Finite  difference  approximations 

5 .  HEIGHT,  TEMPERATURE  AND  VERTICAL  MOTION  PATTERNS  28 

6.  VARIATION  OF  KINETIC  ENERGY  IN  THE  VOLUME  ...  36 

7 .  HORIZONTAL  AND  VERTICAL  ADVECTION  OF  KINETIC  ENERGY  44 

Horizontal  advection  of  kinetic  energy 
Vertical  advection  of  kinetic  energy 


iii 


Page 


8.  GENERATION  OF  KINETIC  ENERGY  AT  THE  BOUNDARIES  52 

Generation  of  kinetic  energy  at  the  vertical 
boundaries 

The  ageostrophic  flux  of  potential  energy 
The  computation  of  and  z' 

Generation  of  kinetic  energy  at  the  boundary 
pressure  surfaces 

Correspondence  between  variations  in  the  (00)) 
process  and  the  u)  patterns 
The  combined  generation  rate  due  to  the  (0vn) 
and  (0to)  processes 

Order  of  magnitude  of  the  generation  rates 

9.  GENERATION  OF  KINETIC  ENERGY  IN  THE  VOLUME  BY 


CONVERSION  FROM  POTENTIAL  ENERGY .  67 

Net  generation  of  kinetic  energy  in  the  volume 

10.  ENERGY  TRANSPORT  AND  GENERATION  PROCESSES  IN 

TERMS  OF  MEAN  AND  EDDY  CONTRIBUTIONS  ....  73 


The  (ku>)  process  -  mean  and  eddy  contributions 
The  (00))  process  -  mean  and  eddy  contributions 
The  (cuT)  process  -  mean  and  eddy  contributions 


11.  KINETIC  ENERGY  BUDGET  OF  THE  VOLUME .  86 

12.  SUMMARY  AND  CONCLUSIONS .  89 

REFERENCES . . .  95 


iv 


LIST  OF  TABLES 

Table  Page 

1.  Kinetic  energy  of  the  volume ,  in  units  of  10  ergs  ...  4l 

24 

2.  Average  values  of  bK/bt,  in  units  of  10  ergs  per  hour  .  43 

3.  Rate  of  change  of  kinetic  energy  of  the  volume  due  to 

advection  of  kinetic  energy, in  units  of  lO2^  ergs  per 

second . . . .  .  45 

24 

4..  Average  values  of  (kvn)  and  (ku>),  in  units  of  10 

ergs  per  hour . .  47 

5.  Rate  of  change  of  kinetic  energy  of  the  volume  due  to 

the  (0vn)  process,  in  units  of  10^1  ergs  per  second  ...  57 

24 

6.  Average  values  of  (0vn),  in  units  of  Hj  ergs  per  hour  .  59 

7.  Rate  of  change  of  kinetic  energy  of  the  volume  due  to 

the  (0iu)  process,  in  units  of  10^4  ergs  per  second  ...  59 

8.  Average  values  of  ($0>),  in  units  of  102^  ergs  per  hour  .  63 

9.  Average  values  of  (0vn)+(0u>),  in  units  of  102^  ergs 

per  hour .  65 

10.  Rate  of  change  of  kinetic  energy  of  the  volume  due  to 
conversion  from  potential  energy,  in  units  of  1021 

ergs  per  second . . .  .  .  68 

/  \  24 

11.  Average  values  of  (cuT),  in  units  of  10  ergs  per  hour  ,  71 

12.  Average  values  of  the  net  generation  rate,  in  units  of 

1024  ergs  per  hour .  72 

13  •  Vertical  transport  of  kinetic  energy  due  to  mean 

motion  and  eddy  components . 78 

14.  Contributions  to  the  (ktu)  process  due  to  mean  motion 

and  eddy  components .  79 


v 


Table  Page 

15.  Vertical  transfer  of  potential  energy  due  to  mean 

motion  and  eddy  components .  8l 

16.  Contributions  to  the  (00))  process  due  to  mean 

motion  and  eddy  components .  82 

17.  Contributions  to  the  (off)  process  due  to  mean 

motion  and  eddy  components .  85 

18.  Average  values  of  terms  in  the  kinetic  energy 

balance  equation  .......  .  87 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1.  Base  map  used  for  analyses  and  grid  points  at  which 

computations  have  been  made . .  .  .  5 

2.  Height,  temperature  and  vertical  motion  at  25  mb, 

50  mb  and  100  mb  at  0300  GOT  on  17  January,  1957  ....  29 

3.  Height,  temperature  and  vertical  motion  at  25  nib; 

50  mb  and  100  mb  at  0300  GCT  on  25  January,  1957  ....  31 

4.  Height,  temperature  and  vertical  motion  at  25  mb, 

50  mb  and  100  mb  at  0300  GCT  on  4  February,  1957  ....  33 

5.  Height,  temperature  and  vertical  nKDtion  at  25  nib, 

50  mb  and  100  mb  at  0300  GCT  on  14  February,  1957  ....  35 

6.  Variation  with  time  of  the  total  kinetic  energy 

in  the  volume . . .  37 

7.  Variation  with  time  of  the  kinetic  energy  in  the 

100-50  mb  layer .  38 

8.  Variation  with  time  of  the  kinetic  energy  in  the 

50-25  mb  layer . .  .  39 

9.  Vertical  transport  rate  of  kinetic  energy . .  50 

10.  Vertical  transfer  rate  of  potential  energy .  62 

11.  Rate  of  conversion  of  kinetic  to  potential  energy  ....  70 


vii 


ABSTRACT 


The  large  and  rapid  decrease  in  kinetic  energy  of  the  lower 
stratosphere  (100-25  mb)  during  the  warming  of  mid- January  to  mid- 
February,  1957  j  h&s  been  investigated  in  terms  of  energy  transport 
and  generation  processes.  These  processes  are  represented  by  the 
corresponding  terms  in  the  mechanical  energy  equation. 

The  evaluation  of  these  terms  was  carried  out  using  the  100- 
rab,  50-mb  and  25-mb  values  of  geopotential,  temperature,  the  geostrophic 
wind  components  and  adiabatic  vertical  motion  at  each  point  of  a  grid 
covering  most  of  the  North  American  continent  and  the  adjacent  North 
Atlantic.  The  investigation  considered  a  limited  volume  which  ex¬ 
changed  energy  with  the  surrounding  atmosphere  through  the  vertical 
boundaries  as  well  as  through  the  upper  and  lower  boundary  pressure 
surface  s . 

Terms  in  the  kinetic  energy  budget  explicitly  involving  exist¬ 
ing  kinetic  energy  were  found  to  be  negligible  in  comparison  with  the 
other  terms.  Conversion  of  energy  within  the  volume  was  from  kinetic 
to  potential  energy  throughout  most  of  the  period. 

Kinetic  energy  was  generated  through  vertical  motions  at  the 
upper  boundary  pressure  surface  and  dissipated  through  vertical  motions 
at  the  lower  boundary  pressure  surface.  The  net  effect  was  a  gain  in 
kinetic  energy  during  the  early  part  of  the  warming  period  and  a  loss 
of  kinetic  energy  during  the  latter  part.  The  contributions  to  kinetic 
energy  changes  on  layers  of  limited  depth  due  to  generation  of  kinetic 
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energy  at  the  upper  and  lower  boundaries  can  be  quite  significant. 

During  the  early  stages  of  the  warming,  the  volume  lost  kinetic 
energy  by  dissipation  through  horizontal  motions  at  the  vertical  bound¬ 
aries.  But  during  the  period  of  pronounced  warming  throughout  the 
volume,  a  large  gain  in  kinetic  energy  was  brought  about  by  generation 
at  the  vertical  boundaries.  This  generation  process  acted  as  a  major 
source  of  kinetic  energy  for  the  volume  and  hence  a  major  sink  of 
kinetic  energy  for  the  lower  stratosphere  at  other  longitudes. 

For  the  volume  and  time  period  chosen  for  study,  large  scale 
sinking  motion  at  the  pressure  surfaces  accounted  for  almost  all  of 
the  energy  transformation  and  vertical  flux  processes. 

The  values  of  the  horizontal  and  vertical  fluxes  of  energy  and 
of  the  total  conversion  rates  were  balanced  against  the  local  change 
in  kinetic  energy  of  the  volume.  The  residual  was  taken  as  an  approxi¬ 
mation  to  the  rate  of  dissipation  of  kinetic  energy  by  small  scale 
motions.  These  rates  appear  reasonable  when  compared  with  those  of  the 
other  processes  in  the  volume. 


1.  INTRODUCTION 


In  recent  years ,  as  a  result  of  increased  frequency  of  observa¬ 
tions  reaching  above  20  km,  several  studies  have  been  made  of  strato¬ 
spheric  circulation  and  in  particular  of  the  warming  phenomena  during 
the  winter  months,  Soherhag  (1952)  first  reported  a  warming  at  strato¬ 
spheric  levels  based  on  ascents  over  Berlin,  Since  then,  other  reports 
of  such  warming  during  successive  winters  have  prompted  the  acceptance 
of  large  temperature  fluctuations  as  a  characteristic  of  wintertime 
stratospheric  circulation.  The  evolution  of  the  warming  as  regards  the 
location  and  intensity  has  varied  from  year  to  year.  But  the  basic 
feature  has  involved  a  fundamental  change  in  the  stratospheric  general 
circulation  from  the  early  winter  type,  consisting  of  a  cold  circum¬ 
polar  cyclonic  vortex,  to  the  spring- summer  type  with  light  winds 
around  a  warm  polar  anticyclone. 

The  stratospheric  warming  in  the  winter  of  1956-1957  was 
particular  significance  not  only  because  of  the  rapidity  with  which 
temperatures  rose  at  levels  above  100  mb  but  also  from  the  point  of 
view  of  its  location  over  North  America.  The  initial  stages  of  the 
warming  (at  least  at  and  below  the  25  mb  level)  took  place  over  south¬ 
eastern  Canada  and  the  subsequent  passage  of  the  warm  center  occurred 
over  an  area  of  relatively  good  data  coverage,  according  to  strato¬ 
spheric  standards.  The  evolution  of  the  temperature  and  wind  fields 
during  this  warming  has  been  described  by  Teweles  (195®)  an^-  by  Craig 
and  Hering  (1959).  Although  such  descriptions  are  of  importance,  an 
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investigation  of  changes  in  the  vertical  motion  patterns  and  associated 
energy  transformations  would  he  of  great  use  in  adding  to  our  knowledge 
of  the  dynamics  of  such  warmings. 

Basic  analyses  of  the  height  and  temperature  fields  at  100  mb, 

50  mb  and  25  mb  during  the  period  1 6  January  to  15  February,  1957  have 
been  carried  out  by  the  stratospheric  circulation  research  group  at  the 
Florida  State  University*  The  period  chosen  includes  a  week  or  so  prior 
to  any  obvious  evidence  of  warming  at  25  mb  and  below,  and  also  extends 
to  a  few  days  after  the  most  marked  temperature  changes  had  occurred. 

The  analysis  area  is  shown  in  Fig.  1.  Both  the  size  of  the  analysis 
area  and  the  upper  limit  to  the  pressure  surface  were  dictated  by  data 
considerations.  Values  of  various  parameters  at  the  three  pressure 
levels  over  a  grid  of  132  points  have  been  made  available  on  punched 
cards.  These  include  values  of  the  geopotential  height,  temperature, 
the  geostrophic  wind  components  and  vertical  motion  as  computed  from 
the  adiabatic  form  of  the  first  law  of  thermodynamics.  Many  of  the 
details  of  the  methods  of  analysis  and  computations  of  vertical  motion 
and  vorticity  values  have  been  described  in  published  reports  (Craig, 
Lateef  and  Mitchem,  196la,  196lb;  Craig  and  Lateef,  1961).  A  summary 
of  the  vertical  motions  has  also  been  presented  by  Craig  and  Lateef 
(1962). 

Studies  of  energy  conversions  at  stratospheric  levels  have  been 
restricted  in  the  past  owing  to  scarcity  of  data  on  vertical  motions  at 
these  levels.  Recent  attempts  to  evaluate  kinetic  energy  generation 
and  dissipation  processes  in  the  lower  stratosphere  were  generally 
confined  to  levels  at  and  below  50  mb.  A  brief  summary  of  previous 


Fig.  1.  Base  map  used  for  analyses  and  grid  points  at  which  computa¬ 
tions  have  been  made.  The  crosses  indicate  intersections  of 
latitude  and  longitude  at  10  degree  intervals,  the  circled 
crosses  indicate  grid  points.  The  perimeter  of  the  grid  is 
shown  in  heavy  lines. 
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studies  on  the  maintenance  of  kinetic  energy  of  stratospheric  motions 
is  given  in  Section  3  of  this  report. 

Although  the  analysis  area  covered  only  a  small  portion  of  the 
Northern  Hemi sphere ,  it  was  felt  that  the  availability  of  data  espe¬ 
cially  on  vertical  motions  up  to  25  mb  and  over  a  sufficiently  dense 
grid  justified  an  investigation  of  the  magnitudes  of  the  terms  that 
enter  into  the  kinetic  energy  budget.  The  computed  values  of  the  vari¬ 
ous  terms  would  indicate  the  relative  importance  of  the  physical  pro¬ 
cesses  contributing  to  the  kinetic  energy  changes  of  the  volume  under 
study.  The  details  of  the  computational  procedures  are  given  in  Sec¬ 
tion  k , 

In  order  to  make  the  results  of  the  computations  more  under¬ 
standable  ,  a  brief  resume  of  the  warming  in  terms  of  the  evolution  of 
the  height  and  temperature  fields  and  of  the  vertical  motion  patterns 
is  presented  in  Section  5. 

The  rest  of  the  report  is  devoted  to  discussion  of  the  results 
of  the  computations.  The  time  variation  of  horizontal  kinetic  energy 
in  the  volume  bounded  by  the  perimeter  of  the  grid  and  the  100 -mb  and 
25-nib  pressure  surfaces  is  presented  in  Section  6.  The  magnitudes  of 
the  various  processes  which  bring  about  the  change  in  kinetic  energy 
are  given  in  the  later  sections. 

In  contrast  to  previous  investigations  of  kinetic  energy  on  a 
hemispheric  basis,  this  investigation  considered  a  limited  volume  which 
exchanged  energy  with  the  surrounding  atmosphere  through  its  vertical 
boundaries  as  well  as  through  the  upper  and  lower  pressure  surfaces. 
Specific  results  concerning  both  the  horizontal  and  vertical,  energy 
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transfers  through  the  boundaries  are  treated  in  Sections  7  and  8, 
while  the  magnitudes  of  the  conversion  processes  between  kinetic  and 
potential  energy  are  discussed  in  Section  9* 

The  limited  longitudinal  extent  of  the  volume  did  not  allow 
evaluation  of  the  energy  transfer  and  conversion  processes  in  terms 
of  motions  on  a  hemispheric  scale.  Analyses  of  the  processes  in  the 
volume  on  the  basis  of  mean  motion  and  departures  from  the  mean,  aver¬ 
aged  both  over  the  area  and  over  suitable  time  intervals  during  the 
period,  was,  however,  undertaken.  The  results  revealed  by  such  analy¬ 
ses  are  included  in  Section  10.  For  the  sake  of  completeness  in  the 
kinetic  energy  budget  of  the  volume,  estimates  of  the  rate  of  frictional 
dissipation  of  kinetic  energy  within  the  volume  were  made  as  described 


in  Section  11. 


2.  THE  MECHANICAL  ENERGY  EQUATION 

Since  on  a  synoptic  scale,  vertical  motion  is  much  smeller  in 
magnitude  than  horizontal  motion,  kinetic  energy  of  horizontal  motion 
can  he  considered  equivalent  to  the  total  kinetic  energy.  The  mathe¬ 
matical  expression  that  represents  the  balance  equation  for  horizontal 
kinetic  energy  can  be  derived  from  the  scalar  product  of  the  horizontal 
wind  and  the  equation  of  motion. 


Derivation 

Using  the  x,  y,  p,  t- coordinate  system,  the  horizontal  equa¬ 
tions  of  motion  may  be  written 


5u 

C>U 

Su 

du 

d0 

-Cr  + 
oc 

U  + 

OX 

v^  + 

‘  S  + 

5v 

dv 

dv 

d0 

*ot 

u^x  + 

V^  + 

= 

fu  +  F 


(1) 

(2) 


where  u  and  v  are  the  horizontal  wind  components  in  the  x  and  y  direc- 
tions  and  u)  =  ™  is  the  individual  pressure  change  representing  vertical 
motion  in  this  coordinate  system  where  pressure  replaces  height  as  the 
vertical  coordinate.  0  is  the  geopotential  of  an  isobaric  surface,  f 
is  the  Coriolis  parameter,  and  Fx  and  F^  are  the  horizontal  components 
of  the  frictional  force  per  unit  mass. 

Multiplying  equation  (l)  by  u  and  equation  (2)  by  v  and  adding 
them,  leads  to  the  following  equation. 
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3k  ,  3k  3k  3k  ,  „ 

St  +  u^+v^+u,^  =  -lu^  +  v 
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2  2 
U  4*  v 

where  k  - - ■? -  ,  is  the  horizontal  kinetic  energy  per  unit  mass. 


Equation  (3)  may  be  re-written  as  follows. 


& 


5u  3v 


Sco 


•  U-f+V 


\ 


+  \u  F  +  v  F 


(4) 


The  last  two  terms  on  the  left  side  of  (4)  cancel  each  other,  since 
by  the  equation  of  continuity,  + 

We  now  select  a  volume  bounded  from  the  surrounding  atmosphere 
by  vertical  boundaries  extending  from  a  constant  pressure  surface,  pQ 
to  another  pressure  surface,  p.  Let  A  denote  the  area  of  intersection 
between  the  vertical  boundaries  and  a  constant  isobaric  surface,  and  L 
be  the  perimeter  of  this  area.  Both  A  and  L  are  considered  independent 
of  height  or  pressure. 

Integrating  (4)  over  the  volume  specified  above,  or  strictly 
speaking  over  the  mass  of  the  atmosphere  between  the  two  constant 
pressure  surfaces  and  the  vertical  boundaries,  leads  to  the  following 
equations. 
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where 


(6) 


is  the  rate  at  which  horizontal  kinetic  energy  of  the  portion  of  the 
atmosphere  being  considered,  is  changing  owing  to  frictional  dissipa¬ 
tion  forces. 

If  we  now  specify  pQ  =  100  mb  and  p  =  25  mb,  we  may  consider 
the  first  term  in  (5)*  namely  -|rr  dA  — ,  to  be  equivalent  to 

p  ^  ^  Pi,  ^  8  p 

•Jr  (°  (  kdA  .  Let  K  denote  the  expression  (  kdA  —  ,  which  is  the 
dt  P  JA  6  Jp  4  g 

total  kinetic  energy  of  the  mass  of  the  portion  of  the  atmosphere  de¬ 
fined  by  A  and  the  100-mb  and  25~mb  preSsurfc  surfaces.  Since  and  p 
have  been  assigned  specific  values,  the  mass  of  the  atmosphere  between 
these  pressure  levels  and  the  vertical  boundaries  remains  constant. 

The  term  c3K/dt,  therefore,  actually  represents  the  time  variation  of 
the  total  kinetic  energy  of  this  constant  mass.  Since  we  ordinarily 
use  data  at  constant  pressure  levels  instead  of  at  fixed  heights  in 
the  atmosphere,  it  has  become  customary  to  speak  of  "volume”  in  the 
x,  y,  p  coordinate  system  corresponding  to  the  geometric  volume  de¬ 
fined  in  the  x,  y,  z  system.  Hence,  the  word  "volume"  wherever  used 
in  this  report  actually  denotes  the  constant  mass  of  the  atmosphere 
between  the  specified  pressure  levels  and  the  vertical  boundaries. 

From  equation  (5),  on  re-arranging  terms,  we  get 
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The  first  term  on  the  right  side  of  (7)  is  the  divergence  of  a 
horizontal  vector  with  components  uk  and  vk.  Therefore,  applying  a  two 
dimensional  form  of  the  divergence  theorem, 


Equation  (7)  can  be  re-written: 


where  the  new  symbols  have  the  following  meanings: 
v^  is  horizontal  wind  component  perpendicular  to  the  vertical 
boundary  of  the  volume,  positive  outward. 

V  is  horizontal  vector  wind. 

V  is  horizontal  del  operator  along  a  constant  pressure  surface. 

The  term  V  •  V0  =  V  .  0V  -  •  V  (10) 


But^  since  V  -  V  * 
we  may  write  (10) 


as 


=  -  a,  where  a  =  —  is  the  specific  volume, 


V  •  ^0  =  V  •  0V  +  (0<u)  + 


mat 


(11) 


Substitution  of  (11)  into  equation  (9)  and  use  of  the  divergence  theo¬ 
rem  finally  yields 
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where  R  is  gas  constant  for  dry  air,  T  is  temperature  and  the  subscript 
p  and  Pq  refer  to  the  upper  and  lower  boundary  values  of  the  correspond¬ 
ing  terms. 

Physical  meaning  of  the  terms 

The  term  dK/dt  on  the  left  side  of  equation  (12)  is  the  time- 
change  of  the  total  horizontal  kinetic  energy  in  the  volume  bounded  by 
the  perimeter  L  and  the  two  specific  pressure  levels  pQ  and  p. 

The  physical  processes  which  bring  about  changes  in  K  are  re¬ 
presented  by  the  terms  on  the  right  side  of  equation  (12)..  Thus,  the 
total  horizontal  kinetic  energy  may  be  changing  as  a  result  of: 

(a)  Advection  of  kinetic  energy  into  the  volume  across  the 
vertical  boundary.  This  is  represented  by  the  term 


(b)  Advection  of  kinetic  energy  into  the  volume  through  the 
boundary  pressure  surfaces.  This  is  represented  by  the  term 
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(c)  Generation  of  kinetic  energy  through  horizontal  motions 
across  the  vertical  boundary.  This  is  represented  by  the  term 

p0  /  ^ 

(0vn)  =  -  |  f  f/L  0vn  dL  I  dp 
4,  N  ' 


(d)  Generation  of  kinetic  energy  through  vortical  motions  at 
the  houndary  pressure  surfaces,,  This  is  represented  hy  the  term 


(e)  Generation  of  kinetic  energy  within  the  volume  hy  conver¬ 
sion  from  potential  energy.  This  is  represented  hy  the  term 


(off)  =  - 


dA  j  dp 

/ 


(f)  Generation  of  kinetic  energy  by  motions  of  smaller  scale 
than  are  included  in  the  analysis.  The  effect  of  such  motions  on  the 
kinetic  energy  of  larger  scales  of  motion  would  ordinarily  amount  to  a 
dissipation  of  the  latter,  and  is  represented  hy  the  term  D. 

Because  of  their  origin  in  the  mathematical  development  of  the 
V  •  V0  term  in  equation  (9);  the  (0vn)  and  (0u))  terms  form  part  of  the 
total  generation  process,  represented  hy  the  integral  of  V  •  ^0.  Over 
the  entire  mass  of  the  atmosphere,  these  two  terms  would  vanish  and  the 
generation  process  would  then  he  given  hy  the  integral  of  ufP,  as  shown 
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in  equation  (12).  For  the  limited  volume,  however,  the  total  generation 
process  may  be  considered  as  the  sum  of  processes  depending  on  the  bound* 
ary  values  of  the  variables,  namely  the  (0v^)  and  (0to)  processes  and 
those  involving  values  within  the  volume  or  the  (off)  process. 

In  the  x,  y,  z,  t  system,  it  can  be  shown  (Starr,  19^8)  that  the 
generation  process  denoted  by  the  integral  of  V  *  Vp  includes  a  (p  v  ) 
term  representing  the  "performance  of  work  by  pressure  forces  at  the 
boundary  in  virtue  of  displacements  due  to  the  horizontal  velocity  com¬ 
ponents."  Inasmuch  as  0  and  u)  in  the  x,  y,  p,  t  system  play  roles  simi¬ 
lar  to  p  and  w  in  the  x,  y,  z,  t  system,  one  can  attribute  similar  inodes 
of  production  of  kinetic  energy  at  the  boundaries  of  the  volume  through 
the  (0vn)  and  (0u>)  processes. 

Under  geostrophic  equilibrium  assumed  for  horizontal  motions 
used  in  the  evaluation  of  a), 

<B=-gpw+pg-||  (13) 

where  w  (=  is  the  vertical  velocity  in  the  x,  y,  z  system.  Equation 
cit 

(13)  indicates  that  the  (0(d)  term  not  only  represents  work  done  on  air 
particles  as  they  move  vertically  through  the  pressure  surfaces  but 
also  involves  release  of  kinetic  energy  through  changes  in  potential 
energy  of  the  atmospheric  layer  due  to  movement  of  the  bottom  or  top  sur¬ 
face  of  the  mass  being  considered. 

To  the  extent  that  the  boundary  generation  processes  are  inter¬ 
linked  with  the  volume  generation  process  (ttfr),  the  (0v^)  and  (0w)  pro¬ 
cesses  are  not  pure  boundary  effects.  They  may,  however,  be  viewed  as 
adjustment  processes  in  response  to  internal  redistribution  of  mass 
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brought  about  by  the  (cut)  process  in  the  volume.  The  specific  nature 
of  generation  of  kinetic  energy  by  the  (0vn)  and  (0U))  terms  as  related 
to  corresponding  potential  energy  changes,  is  not  so  obvious.  However, 
their  inclusion  in  kinetic  energy  generation  rates  in  limited  atmos¬ 
pheric  volumes  would  be  necessary  to  obtain  a  better  estimate  of  the 
total  generation  processes. 

Previous  applications 

In  many  previous  applications  of  equation  (12),  attention  has 
been  focused  on  computation  of  the  (off)  term.  It  is  often  possible  to 
consider  a  volume  that  extends  over  the  hemisphere.  The  (kv^)  and  (0vn) 
terms  are  then  zero,  if  horizontal  transports  across  the  equator  are 
neglected.  If  the  volume  is  based  at  the  earth's  surface  and  extends 
high  enough  so  that  interaction  with  the  overlying  atmosphere  may  be 
neglected,  then  the  (kcu)  and  (0<d)  terms  are  zero, 

Palmen  (i960)  has  discussed  applications  of  equation  (12)  to  a 
volume  extending  from  the  surface  of  the  earth  to  the  top  of  the  atmos¬ 
phere  but  over  a  limited  area.  He  has  called  V  •  V0,  appearing  in 
equation  (9),  the  "work"  term,  representing  work  done  by  the  horizontal 
pressure  forces  inside  the  volume.  For  the  volume  used  by  Palmen,  the 
(0vn)  and  (uoT)  terms  acting  together  constitute  the  generation  process. 

Jensen  (i960)  has  applied  equation  (12)  over  the  Northern  Hemi¬ 
sphere  but  between  the  surface  and  50-mb  pressure  level.  According  to 
him,  the  (0w)  and  (ouT)  terms  together  constitute  the  transformation 
process.  The  transformation  process  now  denoted  by 
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involves  not  only  the  familiar  correlation  between  u)  and  T  but  also  the 
boundary  value  of  the  net  correlation  between  u>  and  0. 

The  dissipation  term  D  is  usually  not  computed  directly,  but 
obtained  as  a  residual  arising  out  of  any  imbalance  among  the  other 
terms.  In  cases  of  computation  of  the  conversion  or  transformation 
process  denoted  by  (off)  alone,  the  value  of  this  term  is  often  taken 
as  an  estimate  of  the  dissipation  rate  of  kinetic  energy  over  a  long 
period  of  time  during  which  dK/5t  is  assumed  to  be  zero. 
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involves  not  only  the  familiar  correlation  between  u>  and  T  but  also  the 
boundary  value  of  the  net  correlation  between  u>  and  0. 

The  dissipation  term  D  is  usually  not  computed  directly,  but 
obtained  as  a  residual  arising  out  of  any  imbalance  among  the  other 
terms.  In  cases  of  computation  of  the  conversion  or  transformation 
process  denoted  by  (tuT)  alone,  the  value  of  this  term  is  often  taken 
as  an  estimate  of  the  dissipation  rate  of  kinetic  energy  over  a  long 
period  of  time  during  which  dK/dt  is  assumed  to  be  zero. 


3.  PREVIOUS  INVESTIGATIONS 


Among  previous  investigations  of  the  energy  budget  of  the  stra¬ 
tosphere,  those  of  White  and  Nolan  (i960),  Wolfe  (i960),  Reed  (1962), 
Jensen  (i960)  and  Barnes  (1962)  are  particularly  relevant  to  the  present 
study. 

White  and  Nolan  (i960)  considered  a  volume  limited  in  the  hori¬ 
zontal  by  latitudes  30N  and  60N  and  longitudes  70W  and  120W  and  bounded 
by  the  200-mb  and  25-mb  pressure  surfaces.  They  studied  the  time  period 
28  January  to  3  February,  1957*  Average  values  of  <jd  and  a  were  obtained 
for  the  layers  200- 50  mb  and  100-25  mb.  Considering  these  values  to  be 
representative  of  tu  and  a  fields  at  the  100  and  50-mb  levels,  they  com¬ 
puted  the  (urr)  term  mentioned  in  Section  2,  namely,  the  conversion  of 
potential  to  kinetic  energy  within  the  volume.  They  found  this  term  to 
be  negative  at  50  mb  that  is,  to  give  conversion  from  kinetic  to  poten¬ 
tial  energy.  The  magnitude  of  the  conversion  process  at  100  mb  was 
nearly  zero.  Although  they  did  not  compute  other  terms  in  the  kinetic 
energy  balance  equation,  they  concluded  that  the  kinetic  energy  of  the 
stratosphere  is  supplied  by  interaction  with  the  adjacent  layers  above 
and  below  the  volume  through  transport  of  existing  kinetic  energy. 

Wolfe  (i960)  investigated  the  rate  of  change  of  kinetic  energy 
of  the  horizontal  wind  and  the  conversion  from  potential  to  kinetic 
energy  at  the  50-mb  level  for  the  warming  period  25  January  to  9  Febru¬ 
ary,  1957*  His  region  comprised  the  latitude  belt  between  35N  and  80N. 
Analysis  was  done  for  only  4  maps  during  the  period.  According  to  him, 
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the  kinetic  energy  of  the  horizontal  wind  showed  a  rapid  and  large 
depletion  due  to  a  large  scale  decrease  of  kinetic  energy  of  the  zonal* 
winds.  The  kinetic  energy  of  the  eddy  flow  showed  a  small,  probably 
insignificant  decrease  during  this  same  period.  The  total  conversion 
of  energy  amounted  to  a  conversion  of  potential  to  kinetic  energy 
although  the  conversion  from  zonal  potential  to  zonal  kinetic  energy 
was  opposite  in  sign.  This  would  contradict  the  results  of  White  and 
Nolan  (i960),  for  the  same  level,  but  for  a  shorter  period  and  over  a 
smaller  region.  An  apparent  discrepancy  arose  between  the  total  rate 
of  change  of  kinetic  energy  and  the  total  rate  of  conversion  of  avail¬ 
able  potential  energy  to  kinetic  energy.  This  was  explained  as  due  to 
treating  the  lamina  as  a  closed  system  and  neglecting  vertical 

flux  of  energy  through  the  lower  boundary  of  the  lamina. 

Reed  (1962)  has  extended  the  previous  investigations  and  com¬ 
puted  the  changes  in  both  the  zonal  and  eddy  components  of  kinetic  and 
potential  energies  as  well  as  corresponding  conversion  rates  for  the 
two  periods  25  January  to  4  February  and  4  February  to  9  February, 

1957*  He  found  a  steady  almost  linear  decline  in  the  zonal  kinetic 
energy  during  the  entire  period.  But  between  25  January  and  b  February, 
while  zonal  kinetic  energy  decreased,  eddy  kinetic  energy  increased  due 
to  conversions  to  and  from  corresponding  potential  energy  sources. 

After  4  February,  however,  these  changes  reversed  their  direction;  a 
decrease  in  eddy  kinetic  energy  occurred  while  there  was  an  apparent 
increase  in  zonal  kinetic  energy  contrary  to  the  decrease  actually  ob¬ 
served.  He  pointed  out  that  conversion  processes  alone  could  not  ac¬ 
count  for  this  discrepancy.  It  appeared  plausible  that  the  loss  of 
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zonal  kinetic  energy  during  the  period  4  February  to  9  February,,  "de¬ 
spite  the  indications  that  other  forms  of  energy  were  being  transformed 
into  this  form  of  energy,  could  be  attributed  to  a  net  vertical  outflow 
of  geopotential  through  the  base  of  the  stratosphere*" 

Jensen  (i960)  carried  out  computations  based  on  data  for  January 
1958  for  about  100  stations  in  the  Northern  Hemisphere  between  20N  and 
80N  in  a  belt  including  all  longitudes*  He  considered  the  volume  from 
the  surface  to  50  mb  and  neglected  any  interactions  with  the  atmosphere 
above.  With  regard  to  results  applicable  to  the  stratosphere,  he  found 
that  the  generation  term  (oil)  was  positive,  although  in  some  layers 
certain  eddy  processes  into  which  this  term  can  be  subdivided,  had  the 
opposite  sign.  More  specifically,  he  found  that  the  stratospheric  layer 
from  200  to  50  mb  was  characterized  by  an  overall  positive  energy  trans¬ 
formation  from  potential  to  kinetic  energy  although  there  was  evidence 
of  a  reverse  transformation  in  the  top  boundary  layer  from  100  to  50  mb. 
The  (0u))  term  was  evaluated  on  the  assumption  that  this  term  was  zero 
at  the  50-mb  surface.  The  contribution  to  the  total  kinetic  energy 
generation  rate  in  the  100- 50  mb  layer,  due  to  this  term  was  found  to 
be  positive.  Neglecting  any  net  transport  in  the  vertical  of  kinetic 
energy  across  the  stratospheric  layers,  dissipation  rates  were  estimated 
as  follows.  It  was  first  assumed  that  the  change  in  kinetic  energy  dur¬ 
ing  the  analysis  period  was  negligible.  Accordingly,  values  correspond¬ 
ing  to  the  (00))  and  (<j)T)  terms  were  introduced  into  the  balance  equation 
(12),  and  the  residual  was  taken  to  represent  the  dissipation  rate. 

Barnes  (1962)  utilized  observations  at  100,  50,  30  and  10-mb 
levels  from  212  stations  in  the  Northern  Hemisphere  for  the  period  July 
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to  December  1957*  He  found  that  in  the  lower  stratosphere,  defined  by 
the  layer  100  to  50  mb,  adiabatic  vertical  motions  resulted  in  an  over¬ 
all  conversion  from  kinetic  to  potential  energy,  in  agreement  with  the 
findings  of  Jensen  (i960).  But  above  to  mb,  the  conversions  seemed  to 
be  in  the  opposite  sense. 

The  (0<d)  term  was  found  by  Barnes  to  be  one  of  the  most  impor¬ 
tant  terms  in  the  kinetic  energy  balance  equation.  Computations  of 
this  term  indicated  that  generation  of  kinetic  energy  at  the  75-mb  and 
to-mb  levels  acted  to  increase  the  kinetic  energy  of  the  layers  above 
these  levels.  Values  at  20  mb  pointed  to  an  increase  in  kinetic  energy 
of  the  layer  below.  But  according  to  him,  these  values  at  20  mb  were 
very  suspect  because  of  the  small  number  of  observations  reaching  10 
mb.  Advection  of  existing  kinetic  energy  was  found  to  be  inadequate  to 
account  for  the  large  increase  in  kinetic  energy  from  October  through 
December  1957* 

More  detailed  studies  are  needed  to  add  to  and  supplement  those 
already  done  for  the  stratospheric  region.  Further  investigations  uti¬ 
lizing  data  above  50  mb  are  required  to  confirm  or  modify  the  conclu¬ 
sions  reached  by  previous  investigators.  While  data  on  a  hemispheric 
basis  would  certainly  be  most  valuable  and  desirable,  there  are  serious 
limitations  on  obtaining  adequate  observations  reaching  even  50  mb. 

The  availability  of  data  reaching  25  mb,  during  the  stratospheric  warm¬ 
ing  of  1957 >  although  for  a  short  period  and  over  a  limited  volume, 
makes  it  possible  to  attempt  a  more  detailed  investigation  of  energy 
transformation  and  flux  processes  in  the  lower  stratosphere.  The  de¬ 
tails  of  the  computational  procedure  used  in  the  evaluation  of  all  the 
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terms  in  the  mechanical  energy  equation  (12)  are  presented  in  the  next 
section. 


k.  COMPUTATIONAL  PROCEDURES 


The  roles  of  the  various  physical  processes  which  affect  the 
balance  of  kinetic  energy  in  a  limited  atmospheric  volume  are  expressed 
by  the  various  terms  in  equation  (12).  If  the  values  of  these  terms 
can  be  separately  determined  from  observational  data,  the  relative  im¬ 
portance  of  the  physical  processes  can  be  ascertained. 

Area  and  line  integrals 

Confutations  were  made  for  an  area  based  on  the  grid  shown  in 
Figure  1.  Lines  were  drawn  midway  between  rows  and  columns  of  grid 
points.  Intersections  of  these  lines  formed  90  squares  with  a  grid- 
point  at  the  center  of  each  square.  Grid  point  values  were  assumed  to 
be  representative  of  the  corresponding  area  increment,  so  that  area  in¬ 
tegrals  were  approximated  by  summation  over  the  grid,  each  value  being 
weighted  by  the  area  involved. 

The  outer  boundaries  of  the  outermost  squares  form  the  perime¬ 
ter  of  the  modified  area.  The  perimeter  consists  of  38  length  incre¬ 
ments,  Under  the  assumption  that  a  value  assigned  to  the  mid-point  of 
a  perimeter  increment  is  representative  of  the  entire  increment,  line 
integrals  were  approximated  by  summation  over  the  perimeter,  each  value 
being  weighted  by  the  length  involved. 

Vertical  integration 

Vertical  integration  through  the  volume  was  approximated  by  the 
trapezoidal  rule: 
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lOOmb 


j  f(p)dp  =  (25mb) 


25  mb 


f(2^)  i  £130} 

2 


+  (50mb) 


f(50)  +  f(lOO) 


=  (I2.5mb)  [f(25)  +  3  f( 50)  +  2  f(lOO)] 


(14) 


Input  data 

Details  of  the  method  of  analysis  and  computational  procedures 
to  obtain  the  values  of  geopotential  height  Z,  temperature  T,  the  geo- 
strophic  wind  components  u  and  v,  and  of  vertical  motion  u)  at  the  lOO¬ 
mb  ,  50-mb,  and  25 -mb  levels  have  already  been  discussed  in  detail  (Craig, 
Lateef,  and  Mitchem,  196la,  1961b;  Craig  and  Lateef,  1961,  1962) ,  Only 
a  brief  summary  of  the  procedures  will,  therefore,  be  given  here. 

The  100-mb  maps  were  analyzed  by  conventional,  subjective  meth¬ 
ods,  including  reference  to  off- time  and  off-level  data  whenever  neces¬ 
sary  and  available.  Careful  attention  was  also  paid  to  the  geostrophic 
wind  relationship.  The  heights  of  the  100-mb  surface  were  considered 
to  be  boundary  values  which,  along  with  the  temperature  distribution  at 
higher  levels  determined  the  topography  of  higher- level  pressure  sur¬ 
faces  according  to  the  hydrostatic  relationship , 

R  ?° 

Z  (x,  y,  p)  =  Z  (x,  y,  pQ)  +  -  T  (x,  y,  p)  d  to  p 

P 

where  Z  is  the  height  of  a  constant  pressure  surface  in  geopotential 
units  and  the  other  symbols  have  the  usual  meanings. 

To  determine  the  mean  temperatures  of  the  100-mb  to  50-mb  and 
50-mb  to  25-mb  layers,  use  was  made  of  time-pressure  cross  sections  pre¬ 
pared  for  all  observing  stations  in  the  analysis  area,  about  130  in  all. 
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The  temperatures  above  100  mb  including  those  pertaining  to  the  50-mb 
and  25-mb  surfaces  were  taken  from  the  analyzed,  smoothed  cross  sections 
and  not  directly  from  the  published  or  original  sources.  The  principal 
advantages  to  this  were,  firstly,  any  diurnal  radiational  errors  were 
smoothed  out  and  secondly,  interpolation  of  missing  data  could  be  made 
with  reasonable  accuracy. 

The  base  map  used  for  all  analyses  was  a  polar  stereographic 
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projection  with  a  scale  of  1:2x10'  at  60N.  The  values  of  height  and 
temperature  were  read  over  the  rectangular  array  of  132  points  (Fig.  l) 
with  a  grid  distance  of  600  km  at  60N  and  of  600  km  (l  +  sin  0)/l.8 66 
at  latitude  0. 

All  computations  have  been  made  for  an  x,  y,  p,  t  coordinate 
system,  with  the  x-y  plane  tangent  to  the  earth's  surface  at  the  loca¬ 
tion  of  the  grid  point  for  which  the  computation  is  made,  and  x  and  y 
positive  in  the  directions  indicated  in  Fig.  1.  Directions  x  and  y  and 
geostrophic  wind  components  u  and  v,  therefore,  point  toward  east  and 
north  only  along  the  central  column  of  grid  points. 

As  with  temperature,  the  errors  in  wind  determination  increase 
with  height  and  are  usually  larger  than  in  the  troposphere,  principally 
because  the  elevation  angle  of  the  sounding  balloon  is  usually  very 
small  by  the  time  the  balloon  reaches  the  tropopause  level.  In  practice, 
considering  the  serious  lack  of  data  above  100  mb  and  the  inaccuracies 
of  the  available  measurements,  it  is  assumed  for  most  purposes  that  geo¬ 
strophic  equilibrium  holds  in  the  stratosphere.  This  assumption  is 
reasonable  in  view  of  the  scale  of  the  systems  involved.  Muench  (1958) 
has  found  the  assumption  to  be  entirely  consistent  with  the  available 
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data. 

The  computations  of  u>  have  been  made  from  the  adiabatic  version 
of  the  first  law  of  thermodynamics  in  the  form, 


<d  = 


ST  &F 
5t  +  u^+  v3£ 


T  dT 

K 


(15) 


where  h  is  the  ratio  of  the  gas  constant  to  the  specific  heat  at  con¬ 
stant  pressure  of  dry  air.  The  accuracy  of  computed  vertical  motions 
in  the  stratosphere,  on  the  assumption  of  adiabatic  temperature  varia¬ 
tions,  has  been  discussed  by  Craig  and  Lateef  (1962). 

With  the  above  approximations  and  with  input  data  consisting  of 
u,  v,  tu,  Z,  T  for  each  grid  point  and  of  appropriate  values  for  area 
and  perimeter  increments,  it  was  possible  to  compute  all  the  terms  in 
equation  (12)  (except,  of  course,  the  dissipation  term). 


Finite  difference  approximations 

The  actual  finite  difference  expressions  used  in  the  computa¬ 
tions  are  as  follows 
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In  the  above  expressions,  AL  and  AA  refer  to  the  perimeter  and  area  in- 
crements  respectively,  v  ,  the  outward  normal  wind  was  evaluated  as 
follows . 

The  average  of  the  v  components  at  grid  points  on  the  first  row 
(Fig.  l)  and  corresponding  points  on  the  second  row  was  taken  as  v^ 
through  the  northern  boundary  of  the  perimeter.  A  similar  procedure 
was  used  for  computing  v^  through  the  southern  boundary.  In  this  case, 
however,  the  v  components  were  used  with  change  of  sign  to  obtain  the 
outward  normal.  Substitution  of  columns  for  rows  and  u 1  s  for  vfs, 
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yielded  vn  at  the  eastern  and  western  boundaries. 

All  computations  were  carried  out  on  the  IBM  709  computer  at 
Florida  State  University. 

Discussions  concerning  the  specific  results  obtained  in  the 
evaluation  of  equations  (l6)  to  (21)  for  twice  daily  maps,  at  0300 
GCT  and  1500  GCT  for  the  period  0300  GOT  of  17  January  1957  to  1500 
GCT  of  Ik  February  1957*  are  contained  in  Sections  6  to  11  of  this 
report. 


5 .  HEIGHT,  TEMPERATURE  AND  VERTICAL  MOTION  PATTERNS 


For  the  purpose  of  making  the  results  of  evaluation  of  the  terms 
described  in  the  previous  section  more  understandable  and  as  a  guide  to 
correlating  variations  in  sign  and  magnitude  of  any  of  the  terms  with 
obvious  changes  in  the  synoptic  pattern,  a  brief  description  of  the 
warming  will  precede  the  presentation  of  the  results.  The  evolution  of 
the  height,  temperature  and  vertical  motion  fields  during  the  course  of 
the  warming  is  described  in  this  section. 

Detailed  presentation  of  the  synoptic  maps,  spaced  1*8  hours 
apart,  of  height,  temperature  and  vertical  motion  are  contained  in  the 
published  reports  already  referred  to  (Craig,  Lateef  and  Mitchem,  196la, 
196lb;  Craig  and  Lateef,  1961) *  Only  a  few  of  these  maps  are,  there¬ 
fore,  shown  in  this  report. 

Figure  2  shows  conditions  on  17  January  1957*  at  the  100-mb, 

50-mb  and  25-nib  surfaces  before  there  was  any  obvious  evidence  of  warm¬ 
ing  at  any  of  these  levels.  A  cold  low  existed  at  all  levels,  somewhere 
near  the  pole,  skirted  by  a  belt  of  strong  westerlies.  A  trough  extended 
southward  from  this  low  across  the  Hudson  Bay  and  over  the  central  United 
States  at  100  mb  and  was  displaced  slightly  to  the  east  at  the  higher 
levels.  The  stratospheric  warm  belt  was  evident  at  100  mb  and  50  mb 
with  relative  temperature  maxima  at  the  trough.  At  25  mb,  the  warm  belt 
was  not  as  evident  as  at  the  lower  levels  due  to  warmer  tenperatures  in 
the  tropical  regions  though  the  temperature  field  showed  two  weak  maxima. 
At  all  levels,  there  was  downward  motion  upstream  from  the  trough  and 
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Fig.  2«  Height,  temperature  and  vertical  motion  patterns  at  0300  GCT 
on  17  January,  1957*  Contours  are  shown  hy  solid  lines  and 
labeled  in  100!s  of  meters  for  100  mb  and  in  10' s  of  meters 
for  50  and  25  mb  (with  the  initial  digit  omitted),  isotherms 
by  dashed  lines  and  labeled  in  °C,  Isopleths  of  u>  are 
labeled  in  units  of  10"^  dyne  cur 2  sec"^.  Cross  hatching  on 
the  zero  isopleth  of  a>  points  toward  downward  nation. 
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upward  motion  downstream  from  the  trough.  Such  a  pattern  of  vertical 
motion  is,  of  course,  what  one  would  expect  if  the  trough  is  a  region 
of  relatively  higher  temperature  and  if  the  time  rate  of  change  of  tem¬ 
perature  is  small  compared  to  the  advection  term  in  the  numerator  of 
equation  (15)- 

During  the  next  few  days  the  trough  at  100  and  50  mb  over  the 
central  United  States,  with  its  associated  temperature  maximum  moved 
slowly  eastward  toward  the  North  Atlantic.  A  new  trough  appeared  at  all 
levels,  over  the  western  United  States  hy  about  21  January.  This  trough 
remained  almost  stationary  until  about  2  February,  when  it  began  to  ret¬ 
rogress  northwestward.  At  100  mb  and  $0  mb,  this  western  trough  was 
warm  compared  to  downstream  temperatures.  In  contrast,  at  25  nib,  tempe¬ 
ratures  increased  in  the  direction  of  flow  along  this  trough  with  no  in¬ 
dication  of  a  relative  maximum  at  the  trough  axis. 

The  pattern  of  vertical  motion  at  100  mb  and  50  mb  during  the 
period  following  17  January  and  until  about  2  February  was  characterized 
by  two  principal  regions  of  ascent  and  descent  in  association  with  the 
two  warm  troughs.  But  the  a)  pattern  at  25  mb  underwent  a  radical  change 
on  and  after  23  January  in  association  with  the  formation  and  intensifi¬ 
cation  of  a  warm  center  at  this  level.  This  warm  center  could  be  located 
off  Nova  Scotia  on  25  January.  Figure  3  shows  conditions  at  all  levels 
on  this  date  and  in  particular  the  occurrence  of  downward  motion  at  25 
mb  over  almost  the  entire  analysis  area.  The  large  scale  downward  mo¬ 
tion  at  25  mb,  initiated  on  23  January  persisted  until  near  the  end  of 
the  analysis  period. 

Changes  in  the  contour  height  field  at  all  levels  and  marked 
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Height,  temperature  and  vertical  motion  patterns  at  0300  GCT 
on  25  January }  1957.  Contours  are  shown  by  solid  lines  and 
labeled  in  100*s  of  meters  for  100  mb  and  in  10* s  of  meters 
for  50  and  25  mb  (with  the  initial  digit  omitted).,  isotherms 
by  dashed  lines  and  labeled  in  °C.  Isopleths  of  u>  are 
labeled  in  units  of  10“1  dyne  cm"2  sec"^.  Cross  hatching  on 
the  zero  isopleth  of  a)  joints  toward  downward  motion. 
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developments  in  the  temperature  and  vertical  motion  fields  at  100  and 
50  mb,  began  to  take  place  after  about  2  February.  The  cold  polar 
vortex  began  to  fill  up  and  move  westward  giving  rise  to  the  establish¬ 
ment  of  east  winds  over  the  northern  portions  of  the  analysis  area.  The 
western  United  States  trough  slowly  retrogressed  northwestward  and  as¬ 
sumed  cold  characteristics  at  all  three  levels,  with  increasing  tempera¬ 
tures  both  upstream  and  downstream. 

The  region  of  coldest  temperature  shifted  from  its  previous  loca¬ 
tion  near  the  arctic  low  to  the  west  or  southwest  of  the  low  center. 
Pronounced  warm  centers  could  be  located  southeast  of  Greenland  at  100 
mb  and  south  of  Iceland  at  50  nib,  while  an  extension  of  the  warming 
region  towards  the  eastern  United  States  was  also  evident.  Figure  4 
illustrates  these  changed  conditions  on  4  February.  By  this  date,  the 
25 -nib  warm  center  which  began  its  intensification  on  23  January  had 
reached  its  maximum  with  a  central  value  of  about  -15C  over  southern 
Greenland.  The  25 -mb  vertical  motions  had  also  reached  their  largest 
magnitudes.  But  radical  changes  had  taken  place  in  the  c 0  fields  at  the 
lower  levels.  As  seen  in  Fig.  4,  the  area  covered  by  downward  motion 
at  these  levels  increased  significantly  and  extended  from  the  cold  west¬ 
ern  trough  to  the  warm  eastern  *  trough.  This  large  scale  motion  per¬ 
sisted  until  near  the  end  of  the  analysis  period. 

During  the  period  4  February  to  the  end  of  the  period,  the  low 
centered  near  Baffin  Island  continued  to  fill  up  rapidly  and  eventually 
showed  up  as  a  weak  cyclonic  circulation  over  the  analysis  area.  The 
axis  of  the  western  United  States  trough  which  was  retrogressing  slowly, 
extended  over  Alaska  towards  the  end  of  the  period.  The  warm  centers 


Fig.  4.  Height,  temperature  and  vertical  motion  patterns  at  0300  GOT 
on  4  February,  1957.  Contours  are  shown  by  solid  lines  and 
labeled  in  100’s  of  meters  for  100  mb  and  in  10' s  of  meters 
for  50  and  25  mb  (with  the  initial  digit  omitted),  isotherms 
by  dashed  lines  and  labeled  in  °C.  Isopleths  of  uj  are 
labeled  in  units  of  10“^  dyne  cm"2  sec"^.  Cross  hatching  on 
the  zero  isopleth  of  cd  points  toward  downward  motion. 
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at  the  lover  levels  also  moved  towards  the  northwest  and  later  toward 
the  west.  By  the  end  of  the  period  (Fig.  5)>  the  north  south  tempera¬ 
ture  gradient  at  all  levels  was  completely  reversed  from  that  prevailing 
in  mid-January.  The  vertical  motion  pattern  at  all  levels  was  marked  by 
two  regions  of  ascent  and  descent ,  generally  over  the  western  and  east¬ 
ern  halves  of  the  continent,  a  reversal  of  the  pattern  on  17  January. 

The  significance  of  the  changes  in  the  vertical  motion  patterns 
that  occurred  at  25  mb  by  about  25  January  and  at  the  lower  levels  after 
about  2  February,  has  been  discussed  by  Craig  and  Lateef  (1962).  They 
have  pointed  out  that  the  important  difference  between  vertical  motions 
prior  to  and  after  about  23  January  at  25  mb  and  before  and  after  2  Feb¬ 
ruary  at  the  lower  levels  lay  in  the  scale  of  downward  motion  at  these 
levels,  rather  than  in  any  large  changes  in  their  magnitude.  The  pre¬ 
ponderance  of  such  motion  at  25  mb  on  25  January  (Fig.  3)  an(3-  its  occur¬ 
rence  at  100  and  50  mb  on  4  February  (Fig.  4)  illustrate  the  large  scale 
of  downward  motion  during  the  stratospheric  warming. 


YlUTrtiL  MOTION 


Fig.  5*  Height,  temperature  and  vertical  motion  patterns  at  0300  GCT 
on  14  February,  1957 «  Contours  are  shown  by  solid  lines  and 
labeled  in  100's  of  meters  for  100  mb  and  in  10' s  of  meters 
for  50  and  25  mb  (with  the  initial  digit  omitted),  isotherms 
by  dashed  lines  and  labeled  in  °C.  Isopleths  of  u>  are 
labeled  in  units  of  10“^.  dyne  cm"2  sec"-*-.  Cross  hatching  on 
the  zero  isopleth  of  u)  points  toward  downward  motion., 


6.  VARIATION  OF  KINETIC  ENERGY  IN  THE  VOLUME 


It  is,  first  of  all,  of  interest  to  consider  how  the  kinetic 
energy  actually  varied  during  the  period.  The  volume  integral  of  ki¬ 
netic  energy  denoted  by  K  on  the  left  side  of  equation  (12)  was  evalu¬ 
ated  twice  daily  for  the  entire  volume  as  well  as  for  the  100- 50  mb  and 
50-25  mb  layers.  Figure  6  shows  the  variation  of  K  with  time  during  the 
period  while  Figures  7  and  8  show  the  variation  of  kinetic  energy  within 
the  layers. 

The  values  of  K  have  been  somewhat  smoothed  by  the  application 
of  an  elementary  binomial  smoothing  function  with  weights  l/4,  l/2,  l/4. 
This  type  of  smoothing  eliminates  fluctuations  in  the  values  due  to  any 
diurnal  variations  and  random  errors.  Following  a  technique  used  by 
Shuman  (1957)  it  can  be  shown  that  the  amplitudes  of  any  variations  of 
wave  components  with  periods  much  longer  than  24  hours  are  not  seriously 
affected  by  this  smoothing  process.  A  smoothed  element  K  defined  as 
above,  is  denoted  by 

^  =  TI  Kt-12  hr  +  2  Kt  +  i  Kt+12  hr  ^ 22  ^ 

With  such  a  definition,  a,  the  factor  by  which  the  amplitudes  of 
the  smoothed  values  of  each  wave  component  are  reduced  is  given  by 

»  ,  1  +  a i  (23) 

where  p.  is  the  wave  number  and  h  is  the  length  of  the  finite  difference 
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TOTAL  KINETIC  ENERGY 
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JANUARY  It67  FEBRUARY 


Fig.  6.  Variation  with  time  of  the  total  kinetic  energy  in  the 
volume ,  bounded  by  the  perimeter  of  the  grid  and  the 
100 -mb  and  25 -mb  pressure  surfaces. 
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Fig.  7»  Variation  with  time  of  the  kinetic  energy  in  the  layer } 
bounded  by  the  perimeter  of  the  grid  and  the  100-rab  and 
50-mb  pressure  surfaces. 


50MB  KINETIC  ENERGY  (l0,4»rg») 
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17  19  21  23  25  27  29  31  2  4  6  8  10  12  14 

JANUARY  1957  FEBRUARY 


Fig.  8.  Variation  with  time  of  the  kinetic  energy  in  the  layer, 
bounded  by  the  perimeter  of  the  grid  and  the  50-rat1  and 
25-mb  pressure  surfaces. 
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increment  in  K.  In  our  case  h  =  12  hours,  since  the  values  of  K  are 
given  12  hours  apart.  According  to  equation  (23),  components  of  wave¬ 
length  10  h,  that  is,  variations  in  K  due  to  fluctuations  with  periods 
of  5  days,  would  be  reduced  only  by  10  per  cent. 

The  time  variations  in  kinetic  energy  exhibited  in  Figs,  6,  7 
and  8  are  almost  identical  in  nature.  Increases  and  decreases  in  kinet¬ 
ic  energy  both  for  the  volume  and  for  the  two  layers  treated  separately, 
occurred  simultaneously.  These  figures  suggest  four  time  periods  into 
which  the  data  might  be  subdivided  for  further  study, 

(a)  From  17  January  to  about  25  January,  K  showed  a  gradual  in¬ 
crease,  with  shorter  period  fluctuations  superimposed.  This  was  the 
period  during  which  a  trough  over  the  central  United  States  moved  east¬ 
wards  to  the  North  Atlantic  while  a  new  trough  appeared  over  the  south¬ 
western  United  States.  The  warming  at  25  mb  had  just  begun  towards  the 
end  of  this  period  but  there  was  no  evidence  of  warming  at  the  lower 
levels. 

(b)  From  26  January  to  about  3  February,  K  showed  a  general  de¬ 
crease.  This  was  the  period  when  the  warming  was  apparent  at  25  mb,  but 

not  yet  throughout  the  entire  volume.  The  cyclonic  center  over  Baffin 
Island  was  slowly  filling  up  at  all  the  three  levels. 

(c)  From  4  February  to  about  8  February,  K  decreased  very  rap¬ 

idly.  This  was  the  period  when  rapid  temperature  increases  were  taking 
place  at  50  mb  and  100  mb.  The  warm  center  at  25  mb  maintained  its  in¬ 
tensity  until  6  February  but  later  started  to  weaken.  The  cold  cyclonic 
vortex  which  remained  nearly  vertical  until  3  February,  began  a  definite 
westward  movement  during  this  period. 
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(d)  From  9  February  to  the  end  of  the  period,  K  continued  to 
decrease,  but  less  rapidly.  This  was  the  period  when  the  rapid  tempera¬ 
ture  changes  had  already  occurred  and  the  north  south  temperature  gradi¬ 
ent  at  all  levels  was  reversed  from  that  prevailing  in  mid-January,  The 
polar  vortex  with  its  belt  of  westerlies  was  replaced  by  a  weak  cyclonic 
circulation  over  the  United  States  with  easterlies  over  the  polar  re¬ 
gions. 

The  values  of  K  for  the  volume  are  given  in  Table  1,  From  these 
values,  the  magnitudes  of  dK/dt  computed  for  the  4  periods  mentioned 
above,  are  given  in  Table  2. 
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Table  1.  Kinetic  energy  of  the  volume,  in  units  of  10  ergs. 


Date 

Time  (GCT) 

K 

17  January 

0300 

127 

1500 

130 

18 

0300 

133 

1500 

140 

19 

0300 

144 

1500 

139 

20 

0300 

135 

1500 

137 

21 

0300 

145 

1500 

147 

22 

0300 

143 

1500 

142 

23 

0300 

143 

1500 

145 

24 

0300 

149 

1500 

157 

25 

0300 

163 

1500 

164 

26 

0300 

164 

1500 

162 

27 

0300 

157 

1500 

156 

Table  1.  Continued 
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Date 

Time  (GCT) 

K 

28  January 

0300 

149 

1500 

138 

29 

0300 

136 

1500 

139 

30 

0300 

i4o 

1500 

144 

31 

0300 

150 

1500 

152 

1  February 

0300 

151 

1500 

143 

2 

0300 

138 

1500 

137 

3 

0300 

137 

1500 

137 

4 

0300 

140 

1500 

136 

5 

0300 

126 

1500 

120 

6 

0300 

115 

1500 

10  4 

7 

0300 

100 

1500 

100 

8 

0300 

91 

1500 

77 

9 

0300 

71 

1500 

69 

10 

0300 

67 

1500 

66 

11 

0300 

66 

1500 

65 

12 

0300 

57 

1500 

52 

13 

0300 

49 

1500 

50 

14 

0300 

51 

1500 
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Table  2.  Average  values  of  dK/dt,  in  units  of  10  ergs  per  hour. 


Period 

dK/dt 

17  January  to  25  January 

0.23 

25  January  to  4  February 

-0.09 

k  February  to  9  February 

-0.57 

9  February  to  Ik  February 

-0.19 

A  comparison  of  these  rates  of  change  in  kinetic  energy  with 

estimates  made  by  Wolfe  (i960)  shows  that  these  values  are  comparable 

to  those  on  a  hemispheric  scale.  Wolfe  estimated  Sk/& t  at  50  nib  to  be 
-2  -1  -1 

about  -3*2  ergs  cm  mb  sec  for  the  period  25  January  to  9  Febru¬ 
ary  1957  over  a  large  portion  of  the  Northern  Hemisphere.  The  corre¬ 
sponding  rate  of  change  of  kinetic  energy  in  the  100- 50  mb  layer  for 

-2  -1  -1 

the  limited  volume  in  the  present  study  is  -2.8  ergs  cm  mb  sec 

No  attempt  has  been  made  in  the  present  study  to  subdivide  the 
kinetic  energy  of  the  volume  into  zonal  and  eddy  components.  With  the 
limited  longitudinal  extent  of  the  volume,  such  an  attempt  is  consid¬ 
ered  neither  useful  nor  meaningful.  It  is  interesting  to  note,  how¬ 
ever,  that  Wolfe’s  computations  showed  that  cHC/dt  was  actually  positive 
in  latitudes  north  of  65N.  The  areas  dominated  by  the  negative  rate  of 
change  of  kinetic  energy  south  of  65N  were  much  larger  than  those  under 
the  influence  of  the  positive  rate  of  change,  with  the  result  that  dK/dt 
assumed  an  overall  negative  value. 


7.  HORIZONTAL  AND  VERTICAL  ADVECTION  OP  KINETIC  ENERGY 


The  horizontal  advection  of  existing  kinetic  energy  across  the 
vertical  boundaries  of  the  volume  is  represented  by  the  (kv^)  berm  of 
Section  2.  Change  in  kinetic  energy  within  the  volume  due  to  vertical 
transport  of  existing  kinetic  energy  through  the  constant  pressure  sur¬ 
faces,  is  given  by  the  (kw)  term. 

Horizontal  advection  of  kinetic  energy 

Parcels  of  air  possessing  kinetic  energy  are  displaced  horizon¬ 
tally  into  or  out  of  the  volume  by  the  wind  components  normal  to  the 
perimeter.  This  leads  to  changes  in  the  kinetic  energy  of  the  volume. 
Computations  of  such  changes,  brought  about  by  horizontal  advection  of 
existing  kinetic  energy,  were  made  according  to  the  finite  difference 
approximation  procedure  shown  in  equation  (17)*  The  method  of  estima¬ 
tion  of  the  values  for  over  each  perimeter  increment  has  already  been 
described  in  Section  4.  Values  of  k  on  the  perimeter  could  not,  however, 
be  computed  because  the  geos trophic  wind  relationship  gave  only  the  u 
components  along  the  first  and  twelfth  rows  of  the  original  grid  of 
132  points  (Pig.  1)  and  only  the  v  components  along  the  first  and  elev¬ 
enth  columns.  In  place  of  the  values  of  k  on  the  perimeter,  the  kinetic 
energies  already  computed  over  the  area  increments  for  use  in  equation 

(16)  were  assigned  to  the  corresponding  perimeter  increments  in  equation 

(17) .  Such  a  procedure  is  satisfactory  enough.  On  the  assumption  of  a 
10  per  cent  error  in  u  and  v  and  a  corresponding  error  in  k  due  to  this 
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procedure ,  it  can  be  shown  that,  in  the  extreme  case  of  v^  being  of  the 
same  sign  all  along  the  perimeter,  the  error  introduced  in  the  evalua¬ 
tion  of  £  (kv  ).  AL.  is  less  than  1  per  cent, 
n'j  J 

The  rates  of  change  of  kinetic  energy  in  the  volume  due  to  hori¬ 
zontal  flux  of  kinetic  energy  across  the  vertical  boundaries,  are  listed 
in  Table  3.  Except  in  the  beginning  of  the  period,  the  horizontal  flux 
consistently  acted  to  decrease  the  kinetic  energy  of  the  volume  until 
about  3  February.  After  this  date  and  until  the  end  of  the  analysis  pe¬ 
riod,  this  process  acted  to  increase  K.  The  rates  of  decrease  in  K  were 
generally  higher  than  the  rates  of  increase. 

Values  of  the  horizontal  advection  term  for  the  100- 50  mb  layer 
and  50-25  mb  layer  at  each  map  time  also  indicated  that  this  term  was 
negative,  that  is,  acted  to  decrease  kinetic  energy,  in  both  the  layers 
up  to  3  February  and  was  positive  thereafter.  These  values  are  not  in¬ 
cluded  in  this  report. 


Table  3* 


Hate  of  change  of  kinetic  energy  of  the  volume  due  to 

21 

advection  of  kinetic  energy,  in  units  of  10  ergs  per  second. 


Bate 

Time  (GCT) 

(*vn) 

(kco) 

17  January 

0300 

0.24 

0.02 

1500 

0.19 

0.03 

18 

0300 

0.04 

0.03 

1500 

-0.16 

0.02 

19 

0300 

-0.20 

0.01 

1500 

-0.13 

0.01 

20 

0300 

-0.16 

0.02 

1500 

-0.24 

-0.01 

21 

0300 

-0.28 

-0.02 

1500 

-0.28 

-0.01 

22 

0300 

-0.25 

0.03 

1500 

-0.23 

0.07 

Table  3.  Continued 
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Date 

Time  (GCT) 

(kvn) 

(koo) 

23  January 

0300 

-0.19 

0.10 

1500 

-0.15 

0.10 

24 

0300 

-0.21 

0.07 

1500 

-0.27 

0.02 

25 

0300 

-0.26 

0.00 

1500 

-0.25 

0.03 

26 

0300 

-0.26 

0.05 

1500 

-0.31 

0.03 

27 

0300 

-0.35 

0.00 

1500 

-0.39 

0.00 

26 

0300 

-0.46 

0.02 

1500 

-0.54 

0.02 

29 

0300 

-0.51 

0.02 

1500 

-o.4i 

0.02 

30 

0300 

-0.34 

0.01 

1500 

-0.29 

0.03 

31 

0300 

-0.26 

0.09 

3.500 

•-O.32 

0.10 

1  February 

0300 

-0.39 

0.07 

1500 

-0.35 

0.05 

2 

0300 

-0.28 

0.01 

1500 

-0.21 

-0.01 

3 

0300 

-0.10 

-0.01 

1500 

0.02 

-0.03 

4 

0300 

0.17 

-0.01 

1500 

0.26 

-0.15 

5 

0300 

0.22 

-0.15 

1500 

0.19 

-0.17 

6 

0300 

0.20 

-0.19 

1500 

0.17 

-0.19 

7 

0300 

0.15 

-0.20 

1500 

0.17 

-0.23 

8 

0300 

0.16 

-0.20 

1500 

0.10 

-0.16 

9 

0300 

0.06 

-0.14 

1500 

0.06 

-0.11 

10 

0300 

0.07 

-0.11 

1500 

0.13 

-O.I3 

11 

0300 

0.16 

-0.14 

1500 

0.16 

-0.12 

12 

0300 

0.15 

-0.10 

1500 

0.10 

-0.08 

13 

0300 

0.05 

-0.08 

1500 

0.04 

-0.09 

14 

0300 

0.04 

-0.09 

1500 

0.03 

-0.08 
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Since  the  u  and  v  components  were  computed  geostrophically,  the 
values  obtained  for  the  horizontal  flux  actually  represent  only  the 
flux  of  geostrophic  kinetic  energy  owing  to  the  geostrophic  wind  com¬ 
ponents.  The  actual  kinetic  energy  may  be  considered  very  nearly  equal 
to  that  computed  with  geostrophic  winds.  But  the  real  flux  may  be  dif¬ 
fer  ent,  as  there  would  be  departures  from  the  actual  wind  both  in  mag¬ 
nitude  and  direction.  Unfortunately;  in  the  absence  of  extensive  wind 
observations;  it  is  impossible  to  know  these  departures.  Estimates  of 
the  magnitudes  of  the  errors  introduced  in  the  computation  of  the  advec- 
tion  term,  using  geostrophic  v^  can  however,  be  made  with  sufficient 
confidence.  Such  estimates  were  not  made  because  this  term  turned  out 
to  be  relatively  unimportant  in  the  kinetic  energy  budget  of  the  volume. 
A  comparison  of  the  magnitudes  of  the  advection  term  with  those  of  the 
generation  terms  is  shown  in  Section  11.  The  extra  effort  in  re-evalu¬ 
ating  the  advection  term  was,  therefore,  considered  unnecessary. 

The  values  of  the  rate  of  change  of  kinetic  energy  in  the  volume 
due  to  horizontal  advection,  for  the  four  time  periods  mentioned  in 
Section  7,  are  contained  in  Table  4. 
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Table  4.  Average  values  of  (kv^)  and  (kw),  in  units  of  10  ergs  per 
hour. 


Period 

(kvn) 

(k(u) 

17  January  to  25  January 

-0.57 

0.12 

25  January  to  4  February 

-1.08 

0.08 

4  February  to  9  February 

0.63 

-0.64 

9  February  to  14  February 

0.33 

-0.38 
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Vertical  advection  of  kinetic  energy 

The  rates  of  change  of  kinetic  energy  in  the  volume  due  to  a 
net  flux  of  kinetic  energy  from  the  adjacent  layers  above  and  below  the 
boundary  pressure  surfaces,  were  computed  according  to  equation  (l8). 
Results  of  the  computations  are  included  in  Table  3. 

The  most  obvious  feature  of  these  results  is  the  occurrence  of 
consistently  negative  values  after  2  February  until  the  end  of  the  pe¬ 
riod.  During  approximately  the  same  period,  the  (kvn)  values  were  con¬ 
sistently  positive  and  of  comparable  magnitude. 

The  values  of  the  vertical  advection  term  were  predominantly 
positive  from  17  January  to  2  February.  The  group  of  small  negative 
values  around  20  January  could,  perhaps,  be  attributed  to  errors  in  data 
particularly  at  the  boundaries  of  the  analysis  area. 

The  values  for  the  rate  of  change  of  kinetic  enegy  in  the  100- 50 
mb  and  50-25  mb  layers  owing  to  vertical  advection  of  existing  kinetic 
energy  above  and  below  these  layers  (not  included  in  this  report)  were 
also  computed  using  equation  (18)  with  the  appropriate  pressure  boundary 
values.  Uhtil  about  31  January,  the  values  for  the  two  layers  treated 
separately,  did  not  show  any  consistent  pattern  as  regards  their  indi¬ 
vidual  contributions  to  the  rate  of  change  of  kinetic  energy  in  the 
volume.  Between  31  January  and  2  February,  the  downward  transport  of 
kinetic  energy  through  the  25-mb  pressure  surface  exceeded  the  amount 
through  the  50~mh  surface  while  the  downward  transport  through  100  mb 
was  less  than  at  50  mb  (Fig.  9).  Thus,  there  was  a  positive  contribu¬ 
tion  to  the  rate  of  change  of  kinetic  energy  in  both  the  layers,  amount¬ 
ing  to  a  substantial  positive  rate  for  the  volume.  The  combined 
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contribution  from  the  vertical  advection  process  through  both  the  layers 
resulted  in  the  group  of  relatively  high  positive  values  in  Table  3  for 
this  period. 

After  about  2  February  and  until  the  end  of  the  period,  the  mag¬ 
nitudes  of  the  downward  transport  of  kinetic  energy  through  the  three 
pressure  surfaces  increased  with  depth  (Fig.  9)«  The  result  was  a  nega¬ 
tive  contribution  to  the  rate  of  change  of  kinetic  energy  in  both  the 
100-50  mb  and  50-25  mb  layers.  For  the  volume,  the  downward  transport 
through  100  mb  exceeded  the  downward  transport  through  25  mb,  so  that 
the  vertical  advection  process  consistently  acted  to  decrease  the  volume 
kinetic  energy.  The  vertical  motion  patterns  at  the  three  pressure  lev¬ 
els,  after  2  February,  show  that  these  negative  values  could  be  largely 
attributed  to  the  radical  changes  that  took  place  in  the  a)  field  at  the 
50  and  100-mb  levels.  As  illustrated  in  Fig.  4,  the  area  of  downward 
motion  at  these  levels  began  to  increase  significantly,  about  this  time. 
The  establishment  of  a  larger  scale  of  downward  motion  at  these  levels, 
rather  than  any  substantial  changes  in  the  magnitudes  of  u),  was  appar¬ 
ently  responsible  for  the  observed  negative  contribution  to  Sk/S t,  due 
to  net  vertical  transport  of  energy  through  the  volume. 

The  values  of  the  rate  of  change  of  volume  kinetic  energy  due  to 
vertical  advection,  for  the  four  chosen  time  periods,  are  listed  in 
Table  4. 

A  comparison  of  the  values  of  (kv^)  and  (kcu)  in  Tables  3  and  k 
indicates  an  apparent  tendency  for  these  two  terms  to  compensate  each 
other.  In  the  presence  of  processes  generating  and  dissipating  kinetic 
energy  in  the  volume,  there  is  no  obvious  reason  to  consider  kinetic 
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JANUARY  1997  FEBRUARY 


Fig.  9.  Vertical  transport  rate  of  kinetic  energy  through  the  100-nib 
(solid  line) ,  50-mb  (dashed  line)  and  25-nib  (dotted  line) 
levels.  Positive  values  signify  downward  transport  of 
kinetic  energy. 
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energy  to  be  locally  conserved.  Any  net  horizontal  advection  of  kinetic 
energy  into  or  out  of  the  volume  through  the  vertical  boundaries,  need 
not,  per  se,  bring  about  a  compensating  vertical  flux  of  existing  kinet¬ 
ic  energy  through  the  horizontal  boundaries.  It  is  interesting  to  note, 
however,  that  Jensen  (i960)  has  also  indicated  the  same  compensating 
tendencies  between  horizontal  convergence  and  vertical  transport  of 
kinetic  energy  in  the  lower  troposphere  over  the  Northern  Hemisphere. 

It  would  appear  therefore,  that  the  apparent  compensation  among  the 
transport  processes  alone,  could  be  a  general  characteristic  of  the 
atmosphere  and  not  a  consequence  of  the  conditions  in  the  limited 
volume  under  study. 


8.  GENERATION  OF  KINETIC  ENERGY  AT  THE  BOUNDARIES 


The  contributions  to  the  rate  of  change  of  kinetic  energy  in  the 
volume  due  to  generating  processes  acting  on  the  boundaries  are  repre¬ 
sented  by  the  terms  (c)  and  (d)  specified  in  Section  2.  Term  (c),  name¬ 
ly,  -  -  (°  ($L  (v  dl)  dp>  gives  the  rate  of  generation  of  kinetic 

g  J  p  L  n 

energy  through  horizontal  motions  across  the  vertibal  boundaries  of  the 

volume.  Term  (d),  namely,  — 

8 

rate  of  generation  of  kinetic  energy  through  vertical  motions  at  the 
boundary  pressure  surfaces. 

Generation  of  kinetic  energy  at  the  vertical  boundaries 

The  rate  of  change  of  kinetic  energy  in  the  volume  due  to  the 

(0v^)  process  was  computed  according  to  the  finite  difference  summation 

procedure  shown  in  equation  ( 1 9 ) *  The  values  of  the  geopotential 

heights  at  the  grid  points  just  inside  the  analysis  area  (Fig.  1)  were 

assigned  to  the  perimeter  increments.  The  error  in  the  evaluation  of 

the  term  £  (Zv  ) .  AL.,  in  equation  (19),  at  any  pressure  level,  by  not 
n.  j  j 

using  Z!s  along  the  perimeter  itself,  amounted  to  less  than  1  per  cent, 
even  under  assumption  of  extreme  situations. 

Initially,  the  geostrophic  v^'s  already  computed  as  described  in 
Section  4,  were  used  for  the  evaluation  of  the  (0vn)  term  for  the  volume, 
for  each  map  time.  This  computation  gave  the  contribution  to  the  (0v^) 
process  from  the  displacements  of  f  at  the  boundary  due  to  the  geo¬ 
strophic  wind  components.  Since  the  generation  process  is  mathematically 


<*») 


“A 


(j  m 


dA 


gives  the 
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equivalent  to  a  flux  or  transfer  rate  of  potential  energy,  the  above 
computation  gave  only  the  gebstrophic  flux  of  potential  energy.  The 
real  flux  might  he  different  to  the  extent  the  actual  wind  departed 
from  the  geostrophic  value.  In  the  absence  of  actual  observed  winds 
at  a  sufficiently  large  number  of  stations  along  the  perimeter,  only 
a  reasonably  close  estimate  could  be  made  of  the  true  flux.  The  meth¬ 
od  of  estimating  the  flux  due  to  the  ageostrophic  component  was  as 
described  below. 


The  ageostrophic  flux  of  potential  energy 

Since  v  can  be  represented  as  v  =  v  -f  v  ,  where  v  and 
n  n  ng  na  ng 

v  are  the  geostrophic  and  ageostrophic  normal  wind  components,  we 
na 

may  write 


(hr  \  =  I  (v  )  +  /  0v  ^ 
/  n  }  ng  /  \  r  na/ 


(24) 


The  ageostrophic  flux  of  potential  energy  is  represented  by  the 


term 


,  .  i  TO  v  *0,  \ 

(  K J-  -  i  j  (#L  dL )  Of  ■  -  j  “•)  <25> 


P  P 

Let  v  %  denote  mean  over  L,  so  that 


-E kl  )  «■  =  E 


r  38  I 

s  (  h  AL  ! 
3=1  J  J ; 


(26) 


and 


0, 


I  0v  \  =  -  L  f  Zv  >  dp 
y  na  J  j  na  * 


(27) 


The  evaluation  of  Zv^a  required  knowledge  of  v^a  over  each  peri¬ 


meter  increment.  This  was  not  available  in  the  absence  of  actual  wind 


observations  along  the  perimeter.  The  following  procedure  was  therefore 


adopted  to  obtain  a  close  estimate  of  (0vna)  for  the  volume. 

The  averaging  of  the  values  of  the  product  of  Z  and  along 
the  perimeter  can  also  be  expressed  as 


Zv  =  Z 
na 


v  + 
na 


Z’v! 


na 


(28) 


where  the  primed  quantities  denote  departures  of  Z  and  v^  over  each 
perimeter  increment,  from  their  peripheral  mean  values. 

- - 1  V”-'\ 

In  what  follows,  the  method  of  obtaining  the  mean  term  Z  v^a  is 

described.  It  will  be  shown  that  the  term  Z’v'  could  be  neglected  as 

na 

being  at  least  an  order  of  magnitude  smaller  than  the  mean  term. 


The  computation  of  v  ^  and  Z 
na 


From  the  relations 


Lv  =  L  v  +Lv 
n  ng  na 


?0  ,  ?0  x  /  r 

I  1?*  VdAdp  =  -  j  ^  |  j 
p  A  p  '  A 


CO  dA  I  dp 


and 


j  V  *  V  dA  =  0T  v  dL  =  L 
J  -  j  L  n  n 


(29) 

(30) 


(31) 


it  follows  that  v^  can  be  computed  from  v^  given  by  equations  (30)  and 

(31)  and  from  v  ^  obtained  through  the  geostrophic  u  and  v  components 
ng 

across  the  perimeter  increments. 

However,  equation  (30)  implies  that  the  values  of  v^  would  essen¬ 
tially  apply  to  the  mean  pressure  height  of  the  layer  through  which  the 
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integration  is  carried  out.  Hence  for  the  100-50  mb  layer  and  v^ 

for  the  50-25  mb  layer  were  computed  separately  and  these  values  of 

v^  were  considered  to  apply  at  75-^b  and  37.5-mb  levels  respectively. 

75  mb  was  taken  to  be  the  mean  of  at  100  mb  and  50  mb  and 

v  at  37*5  mb  to  be  the  mean  of  values  at  50  and  25  mb.  Using  the 
ng  w 

values  of  ^  and  vj  at  75  and  37-5  mb,  (v^),^  and  were 

computed  for  each  map  time. 

The  next  step  was  to  obtain  the  geopotential  height  Z  at  75  mb 

and  37.5  mb.  On  the  assumption  of  a  logarithmic  relationship  between 

height  and  pressure,  the  following  approximations  were  used  to  get  Z^^ 

and  Z__  _  ,  . 

37 • 5mb 


7  -  7  ^  +  in  10  ^  xn  7*5  i 

75mb  ~  lOOmb  ^  in  10  -  in  5  !  50mb 


"'lOOmb  : 
1 


=  O.585  Z,—  .  +  0.415  Zc  '  . 

lOOmb  y  50mb 


(32) 


Similarly, 


*  0-585  Vinb  *  0'ld5  5«a> 


(33) 


Having  obtained  Z  and  v  at  75  rob  and  37*5  mb,  the  contribution 

na 

to  the  (0vn)  process  of  the  term  L  Z  v  could  be  computed. 


For  the  100- 50  mb  layer,  this  term  amounted  to  50  mb 


(zr'l 

\  na  1 


na  ^75mb  for  the  5^-25  mb  layer,  it  amounted  to  25  mb 
L  ^Z  (_rWK.  The  sum  of  these  two  terms  gave  an  estimate  of  the 


na 


-j  37  •  5®b  * 


contribution  to  the  (0vn)  process  due  to  the  ageostrophic  wind  compo¬ 


nents. 


Though  the  term  Z’v'  of  (28)  could  not  be  evaluated  in  the 
na 

absence  of  actual  wind  observations,  it  is  possible  to  estimate  the 
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relative  order  of  magnitude  of  this  term  as  compared  to  Zv^. 

> — , 

As  part  of  the  computations  of  Z ,  the  values  of  Z*  over  each 

perimeter  increment  were  also  obtained.  These  values  indicated  that 

the  maximum  value  of  Z1  at  any  level  did  not  exceed  500  meters ,  during 

the  period.  Reference  to  the  order  of  magnitude  of  the  values  of 

L  Zv  term  at  any  level,  showed  that  these  values  were  of  the  order 
ng  J  ’ 

l8  O'  i  *  ■ ...  — — 

of  10  cnr  sec  .  Hence,  in  order  that  Zfv'  should  make  a  sizeable 

’  na 

contribution  to  the  (0v^)  process,  the  integrated  value,  L  Z'v1^, 

should  be  of  the  same  order  of  magnitude.  Knowing  the  maximum  value  of 

Zf,  the  magnitude  of  v*  corresponding  to  a  value  for  L  Z*vf  of  say, 

na  na 

17  3  -1 

10  cm  sec  ,  can  be  ascertained. 


38/ 


Let  AL  =  5^0  km,,  an  average  value.  Then  L  Z'v' 


na 


E  f Z 1  v '  \  .  AL.  =  1017  cm3  sec"1,  or 

j=l'  na.)  J  J 


38  (5  x  10^  cm)  (v1  )  (5.L  x  107  cm)  =  lO'*’7  cm3  sec”\ 

na 

Thus  vf  should  be  at  least  10  meters  per  second,  over  every 
na 

perimeter  increment,  each  one  of  them  having  Z1  =  500  meters.  It  is 


difficult  to  visualize  a  departure  of  v'  of  this  magnitude  acting  in 

na 

the  same  direction  all  around  the  perimeter  at  all  levels,  together 
with  the  unrealistic  extreme  departure  Zf  of  500  meters  over  every  peri¬ 
meter  increment.  The  contribution  of  Z-'v'  to  the  (0vn)  process  may, 
therefore,  be  neglected. 

Table  5  gives  the  values  of  the  (0vng)  term,  the  (0vnQ)  term  and 
the  sum  of  these  two  terms  representing  the  rate  of  change  of  kinetic 
energy  in  the  volume  due  to  the  (0vn)  process. 
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Table  5.  Rate  of  change  of  kinetic  energy  of  the  volume  due  to  the 

21 

(0vn)  process,  in  units  of  10c  ergs  per  second. 


Date 

Time  (GCT) 

(0v  ) 
r  ng' 

(0v  ) 

na' 

(0vn) 

17  January  0300 

-16.85 

7.96 

-8.89 

1500 

-10.53 

3.98 

-6.55 

18 

0300 

-5.69 

0.1 

-5-59 

1500 

-1.45 

-1.71 

-3.16 

19 

0300 

3.51 

-3.63 

-0.12 

1500 

7.61 

-7.03 

0.58 

20 

0300 

9-57 

-8.35 

1.22 

1500 

11.91 

-8.25 

3*66 

21 

0300 

14.54 

-11.15 

3-39 

1500 

16.57 

-16,20 

0.37 

22 

0300 

16.91 

-18.48 

-1.57 

1500 

13.36 

-16.43 

-3.07 

23 

0300 

8.45 

-13.29 

-4.84 

1500 

6.69 

-II.87 

-5.18 

24 

0300 

8.45 

-11.06 

-2.61 

1500 

9.94 

-8.26 

1.68 

25 

0300 

7.69 

-4.71 

2.98 

1500 

4.47 

-4.20 

0.27 

26 

0300 

3-77 

-6.53 

-2.76 

1500 

2.62 

-6.90 

-4.28 

27 

0300 

-3.09 

-2.22 

-5-31 

1500 

-9.96 

2.31 

-7.65 

28 

0300 

-10.15 

0.51 

-9.64 

1500 

-6.24 

-5-53 

-11.77 

29 

0300 

-3.34 

-4.85 

-8.19 

1500 

-1.94 

-1.96 

-3.90 

30 

0300 

-0.99 

-3.72 

-4.71 

1500 

-1.89 

-5.42 

-7.31 

31 

0300 

-6.26 

-2.74 

-9.00 

1500 

-6.56 

-2.57 

-9.13 

1  February  0300 

-2.58 

-5.77 

-8.35 

1500 

-O.31 

-6.90 

-7.21 

2 

0300 

2.44 

-5.23 

-2.79 

1500 

8.37 

-5-75 

2.62 

3 

0300 

17.09 

-11.55 

5.54 

1500 

27.38 

-17.84 

9-5^ 

4 

0300 

35-49 

-18.60 

16.89 

1500 

38.80 

-15.63 

23.17 

5 

0300 

40.75 

-14.45 

26.30 

1500 

^3-73 

-14.14 

29.59 

6 

0300 

45.43 

-11.88 

33.55 

1500 

44.89 

-9.89 

35.00 

7 

0300 

45.29 

-10.26 

35.03 

1500 

45.11 

-8.11 

37.00 
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Table  5.  Continued 


Date 

Time  (GCT) 

(0V  ) 

ng 

(SZSv  ) 
na7 

8  February 

0300 

40.59 

-1.26 

39  •  43 

1500 

34.90 

3-92 

38.82 

9 

0300 

29.26 

5.19 

34.45 

1500 

25.67 

3.52 

29.19 

10 

0300 

27.22 

0.77 

27.99 

1500 

29.76 

2.16 

31.92 

11 

0300 

30.09 

5.61 

35-70 

1500 

28.15 

6.95 

35-10 

12 

0300 

23.73 

8.38 

32.11 

1500 

20.49 

10.95 

31.44 

13 

0300 

18.89 

14.46 

33-35 

1500 

17.72 

15.99 

33.71 

14 

0300 

17.85 

14.43 

32.28 

1500 

18.07 

11.88 

29.95 

It  is  evident  from  these  values  that  the  (0v^)  term  computed 
with  geostrophic  winds  alone,  could  be  very  much  in  error.  Further, 
since  the  magnitude  of  this  term  is  largely  dependent  on  0,  any  errors 
in  v^  due  to  the  geostrophic  assumption  would  be  magnified  and  could 
lead  to  spurious  dissipation  rate  estimates  in  the  final  energy  balance, 
So  far  as  can  be  ascertained,  no  previous  attempts  at  evaluating  this 
term  for  a  limited  volume  have  been  made.  On  the  other  hand,  in  pre¬ 
vious  investigations  (Jensen,  i960;  Barnes,  1962)  of  the  kinetic  energy 
budget  of  the  stratosphere  over  the  Northern  Hemisphere,  either  this 
term  was  neglected  or  only  a  partial  evaluation  of  its  magnitude  at  the 
equator  was  carried  out. 

The  rates  of  change  of  kinetic  energy  in  the  volume  during 
the  four  sub-periods  due  to  the  (0vn)  process  are  shown  in  Table  6. 
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Table  6.  Average  values  of  (0vn) ,  in  units  of  102^  ergs  per  hour. 


Period 

(0Vn) 

17  January  to  25  January 

-5.91 

25  January  to  4  February 

-Ilk  77 

4  February  to  9  February 

116.86 

9  February  to  14  February 

116.55 

Generation  of  Kinetic  energy  at  the  boundary  pressure  surfaces 

The  rates  of  change  of  kinetic  energy  in  the  volume  as  well  as 
in  the  100 -50  mb  and  50-25  mb  layers  due  to  the  (0u))  process  were  com¬ 
puted  according  to  the  finite  difference  summation  procedure  shown  in 
equation  (20).  The  results  of  the  computations  are  shown  in  Table  7. 

Table  7  •  Hate  of  change  of  kinetic  energy  of  the  volume  due  to  the 
(0uj)  process,  in  units  of  10^  ergs  per  second. 


Date 

Time  (GCT) 

(0<o) 

100 -50  mb 

(00)) 

50-25  mb 

(00)) 

Volume 

17  January 

0300 

-2.14 

5-57 

3.44 

1500 

0.87 

4.48 

5.35 

18 

0300 

1.10 

3.49 

4.59 

1500 

-1.88 

3.38 

1.51 

19 

0300 

-5.16 

4.08 

1.08 

1500 

-3.36 

4.95 

1.58 

20 

0300 

-O.19 

4.31 

4.12 

1500 

-2.82 

4.05 

1.23 

21 

0300 

-3.02 

5.11 

2.09 

1500 

3.51 

4.18 

7.68 

22 

0300 

5.90 

2.48 

8.38 

1500 

7.68 

2.83 

10.51 

6o 


Table  7*  Continued 


Date 

Time  (GCT) 

(0m) 

100-50  mb 

(M 

50-25  mb 

(00)) 

Volume 

23  January  0300 

1.09 

3-  -5 

14.11 

1500 

1.12 

31^8 

14.86 

24 

0300 

8.72 

4*'  3 

13.45 

1500 

3.98 

4.i  4 

8.22 

25 

0300 

3.30 

4.C  5 

7-35 

1500 

7.17 

4.9  *■ 

12.11 

26 

0300 

10.62 

4.9  • 

15.57 

1500 

13.43 

4.0  ; 

17.45 

27 

0300 

12.55 

4.9' 

17.52 

1500 

11.88 

7-5S 

19.47 

28 

0300 

15.65 

7-7^ 

23.40 

1500 

17.67 

7-03 

24.71 

29 

0300 

13.93 

9.12 

23.05 

1500 

8.73 

12.56 

21.29 

30 

0300 

5.89 

13.54 

19.43 

1500 

7-35 

11.97 

19.33 

31 

0300 

12.03 

11.01 

23.04 

1500 

15.38 

9.53 

24.90 

1  February  0300 

18.53 

7-57 

26.10 

1500 

17.16 

7-79 

24.96 

2 

0300 

17.28 

6.77 

24.05 

1500 

8.78 

2.63 

11.41 

3 

0300 

10.76 

0.81 

11.57 

1500 

8.36 

2.47 

10.83 

4 

0300 

1.98 

2.01 

3-99 

1500 

-1.50 

-1.22 

-2.72 

5 

0300 

-0.64 

-4.61 

-5.25 

1500 

-1.36 

-6.99 

-8.35 

6 

0300 

-3.60 

-7.49 

-11.10 

1500 

-3-44 

-8.87 

-12.31 

7 

0300 

-1.71 

-12.52 

-14.23 

1500 

-1.95 

-15.90 

-17.85 

8 

0300 

-1.16 

-17.63 

-18.79 

1500 

-1.91 

-18.86 

-20.77 

9 

0300 

-5.30 

-13.27 

-18.58 

1500 

-0.99 

-14.44 

-15.43 

10 

0300 

-0.95 

-11.43 

-12.38 

1500 

-7-03 

-11.98 

-19.01 

11 

0300 

-11.04 

-13.22 

-24.26 

1500 

-11.60 

-11.95 

-23.54 

12 

0300 

-9.79 

-10. 31 

-20.09 

1500 

-9.24 

-10.67 

-19.91 

13 

0300 

-14.02 

-10.42 

-24. 44 

1500 

-16.65 

-8.47 

-25.12 

14 

0300 

-14.37 

-8.20 

-22.58 

1500 

-12.38 

-10.61 

-22.99 
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The  most  striking  feature  of  these  results  is  the  occurrence 
of  consistently  positive  values  for  the  volume  from  the  "beginning  of 
the  period  until  about  4  February  and  consistently  negative  values 
thereafter.  In  terms  of  rate  of  change  of  kinetic  energy  in  the  volume, 
this  implied  that  the  volume  gained  kinetic  energy  during  the  period 
17  January  to  4  February  and  lost  kinetic  energy  after  4  February, 
through  net  generation  at  the  upper  and  lower  boundary  surfaces.  It  is 
also  of  interest  to  note  that,  except  in  the  beginning  of  the  period, 
the  volume  always  gained  kinetic  energy  through  generation  at  the  25 -mb 
surface  and  lost  kinetic  energy  through  dissipation  at  the  100-mb  sur¬ 
face.  This  can  be  seen  from  Fig.  10  in  which  the  values  of  £  ( Zco) ^  AA^, 
at  each  level,  are  plotted  for  each  map  time. 

On  the  first  two  days,  the  generation  process,  mathematically 
equivalent  to  the  vertical  transfer  of  potential  energy,  showed  upward 
transfer  through  the  three  pressure  surfaces.  This  implied  gain  of 
kinetic  energy  by  generation  through  net  upward  motion  at  100  mb  and 
loss  by  dissipation  through  net  upward  motion  at  25  mb.  The  gain  was 
greater  than  the  loss,  so  that  the  (0u>)  term  was  positive. 

From  19  January  until  the  end  of  the  period,  the  values  of 
E  (Zu))^  AA^  at  each  level  were  consistently  positive,  implying  a  down¬ 
ward  transfer  of  potential  energy.  During  the  period  19  January  to  4 
February,  the  gain  of  kinetic  energy  due  to  generation  through  net  down¬ 
ward  motion  at  25  mb  exceeded  the  loss  due  to  dissipation  through  net 
downward  motion  at  100  mb.  After  4  February,  the  loss  through  dissipa¬ 
tion  at  100  mb  exceeded  the  gain  through  generation  at  25  mb. 

The  rate  of  generation  of  kinetic  energy  in  the  5°-25  mb  layer 


POTENTIAL  ENERGY  TRANSFER  RATE  ( IO*Wg  we"') 
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JANUARY  1957  FEBRUARY 


,  Vertical  transfer  rate  of  potential  energy  through  the 
100 -mb  (solid  line)  50-mb  (dashed  line)  and  25-mb 
(dotted  line)  surfaces.  Positive  values  denote  down¬ 
ward  transfer  of  energy. 


Fig.  10 
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was  positive  before  4  February  and  negative  thereafter.  The  ( jZ5u>)  pro¬ 
cess  in  the  100- 50  mb  layer  maintained  the  same  sign  as  in  the  higher 
layer  from  about  21  January  till  the  rest  of  the  period.  The  reversal 
in  sign  that  occurred  on  4  February  was  simultaneous  in  both  the  layers. 

The  rates  of  change  of  kinetic  energy  in  the  volume  during  the 
four  subperiods  due  to  the  (0tu)  process,  are  shown  in  Table  8. 

Table  8.  Average  values  of  (0U)),  in  units  of  10^  ergs  per  hour. 


Period 

(00)) 

17  January  to  25  January 

24.86 

25  January  to  4  February 

67.74 

4  February  to  9  February 

-38.67 

9  February  to  14  February 

-73-82 

Correspondence  between  variations  in  the  (0to) 

process  and  the  u>  patterns 

It  is  of  interest  to  view  the  variations  in  both  magnitude  and 
sign  of  the  (0u))  process  as  related  to  the  radical  changes  that  occurred 
in  the  cu  fields  at  the  three  levels  during  the  course  of  the  warming. 


Except  during  the  first  few  days,  the  positive  rate  of  genera¬ 
tion  of  kinetic  energy  in  the  volume  as  well  as  in  the  two  layers  before 
4  February,  was  brought  about  by  a  larger  generation  rate  or  synonymously 
speaking,,  a  larger  downward  transfer  of  potential  energy  through  the 
upper  boundary  (25  mb)  surface  as  compared  to  the  downward  transfer 
through  the  lower  boundary  surface  (100  mb).  The  cu  maps  for  25  January 
indicate  that  the  large-scale  downward  motion  at  25  mb,  contributed  to 
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such  a  large  downward  transfer  at  this  level.  *  The  occurrence  of  large 
scale  downward  motion  at  25  mb  from  23  January  to  about  2  February  with 
no  marked  changes  in  the  a)  patterns  at  the  lover  levels,  was  associated 
with  the  positive  values  of  the  ($u>)  term  or  relatively  high  rates  of 
generation  of  kinetic  energy  during  this  period. 

After  2  February,  radical  changes  took  place  at  the  50“mt)  and 
100-mb  levels.  As  described  in  Section  5/  there  was  a  tendency  for 
downward  motion  at  these  levels  to  occupy  increasingly  larger  areas 
(Fig.,  4),  Associated  with  this,  the  downward  transfer  of  energy  through 
the  100-mb  surface  began  to  increase  in  magnitude.  By  about  4  February, 
the  downward  transfer  through  100  mb  exceeded  that  through  the  higher 
levels,  resulting  in  large  negative  values  of  the  (0co)  term  or  dissipa¬ 
tion  of  kinetic  energy  in  the  volume. 

The  changes  in  the  a)  fields  and  the  associated  variations  in 
magnitude  and  sign  of  the  energy  generation  rate  in  the  volume  through 
generation  at  the  boundaries,  can  be  briefly  summarized  as  follows. 
Generation  of  kinetic  energy  at  the  25-mb  surface  in  association  with  a 
general  sinking  motion  at  this  level  acted  to  increase  the  kinetic 
energy  of  the  underlying  atmospheric  volume.  The  establishment  of 
large  scale  sinking  motion  at  the  lower  pressure  levels  acted  to  de¬ 
crease  the  kinetic  energy  of  the  overlying  layers  through  dissipation 
at  the  lower  boundary  surfaces. 

The  combined  generation  rate  due  to  the  (0vn)  ar*d  (0U))  processes 

For  convenient  reference,  the  values  of  the  (0v^)  and  (0u>)  terms 
and  the  total  generation  rate  they  represented,  are  given  below. 
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pk 

Table  9.  Average  values  of  (0vn)+(0u)),  in  units  of  10  ergs  per  hour. 


Period 

(0vn) 

CfM 

(0vn)+($M 

17  January  to  25  January 

-5.91 

24.86 

18.95 

25  January  to  4  February 

-14.77 

67.74 

52.97 

4  February  to  9  February 

116.86 

-38.67 

78.19 

9  February  to  14  February 

116.55 

-73.82 

42.73 

As  in  the  case  of  the  advection  processes  described  in  Section 
1,  the  (^v^)  and  (0(u)  processes  also  acted  in  opposite  sense  or,  par¬ 
tially  compensated  each  other  to  give  a  net  positive  kinetic  energy  gen¬ 
eration  rate  for  the  volume.  In  terms  of  the  dynamics  of  the  processes, 
the  loss  of  kinetic  energy  through  dissipation  at  the  vertical  boundaries 
and  at  the  100-mb  surface,  was  compensated  by  a  gain  through  generation 
at  the  25-mb  surface.  Later,  a  gain  in  kinetic  energy  through  generation 
at  the  vertical  boundaries  and  at  the  25-mb  surface,  was  reduced  by  dis¬ 
sipation  at  the  100-mb  surface.  If  we  view  the  processes  as  equivalent 
to  flux  of  potential  energy,  horizontal  divergence  or  convergence  of  the 
potential  energy  flux  produced  a  corresponding  vertical  convergence  or 
divergence  of  the  flux  through  the  constant  pressure  surfaces. 

Order  of  magnitude  of  the  generation  rates 

The  values  obtained  for  the  (0vn)  and  (0co)  processes  are  an  order 
of  magnitude  larger  than  the  corresponding  values  for  (kv^)  and  (kuu)  pro¬ 
cesses  discussed  in  the  previous  section.  This  is  in  agreement  with 
what  one  would  expect  by  a  scale  factor  analysis.  Denoting  scale 
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equivalence  by  the  symbol  [  ],  we  may  write  V  =  j  ^ 

is  the  wind  velocity  and  L  is  a  characteristic  length. 


\L 

fL 


,  where  V 


kinetic  energy 

The  ratio  — - — varies  as 
potential  energy 


Relative  vorticity 


i  ji.  jl 

^  I  f  f  L2 


i  = 


'  V 

= 

'JL' 

_LJ 

.fi2. 

Thus, 


4 


UM  = 


[1/10] 


(34) 


(35) 

(36) 


As  mentioned  in  Section  7 >  the  (kv^)  term  representing  horizontal 
flux  of  kinetic  energy  was  evaluated  using  the  geostrophic  u  and  v  com¬ 
ponents  across  the  perimeter  increments.  In  view  of  the  near  geostrophic 
equilibrium  observed  at  [stratospheric  levels  and  the  fact  that  the  (kv  ) 
term  was  relatively  unimportant  in  the  kinetic  energy  balance,  the  evalu¬ 
ation  of  (kv  )  term  was  not  considered  essential.  But  the  inclusion  of 
na 

the  (0vna)  process  in  the  evaluation  of  the  (0v^)  term,  was  necessary  in 
order  to  obtain  a  closer  estimate  of  the  actual  horizontal  flux  of  poten¬ 
tial  energy.  It  is  easily  seen  that  the  integral  of  0v^  around  the 
closed  perimeter  L  would  be  zero  but  for  the  variations  of  f  along  L. 

The  (0v  )  term,  however,  depended  heavily  on  the  magnitudes  of  0.  Any 

na 

errors  in  v^  due  to  the  ageostrophic  component,  would,  therefore  be 
greatly  magnified  in  the  summation  of  the  products  of  0  and  v  and  could 
lead  to  spurious  values  for  the  (0v^)  term. 


9.  GENERATION  OF  KINETIC  ENERGY  WITHIN  THE  VOLUME  BY  CONVERSION  FROM 
POTENTIAL  ENERGY 

The  rate  of  conversion  of  available  potential  energy  into  the 
kinetic  energy  of  the  horizontal  wind  is  given  by  the  ((UP)  term  of  equa¬ 
tion  (12).  The  finite  difference  approximation  used  for  the  evaluation 
of  this  term  is  shown  in  equation  (21)*  That  this  term  actually  repre¬ 
sents  a  conversion  of  potential  and  internal  energy  to  kinetic  energy 
has  been  demostrated  by  previous  investigators  (e.g.,  White  and  Saltz- 
man,  1956).  The  physical  mechanism  for  such  conversion  is  the  sinking 
of  cold  air  mass,  and  the  rising  of  warm  air  mass.  In  what  follows, 
any  statement  about  conversion  of  potential  energy  implies  conversion 
of  the  sum  of  potential  and  internal  energy  of  the  volume. 

The  rates  of  generation  of  kinetic  energy  in  the  volume  due  to 
the  conversion  process,  (<0T),  are  given  in  Table  10.  Except  in  the  be¬ 
ginning  of  the  period,  this  process  consistently  represented  dissipa¬ 
tion  of  kinetic  energy  in  the  volume  through  conversion  to  potential 
energy.  The  values  for  the  conversion  term  within  the  100- 50  mb  and 
50-25  mb  layers  treated  separately,  are  plotted  in  Fig.  11.  During  the 
first  two  days  of  the  period,  kinetic  energy  in  both  the  layers  was  in¬ 
creasing  due  to  conversion  from  potential  energy.  From  19  January  to 
about  24  January,  the  kinetic  energy  dissipation  rate  due  to  a  reverse 
transformation  from  kinetic  to  potential  energy  was,  in  general,  of  the 
same  magnitude  in  both  the  layers.  The  dissipation  rate  due  to  the  con¬ 
version  process  was  higher  in  the  50-25  mb  layer  than  in  the  100-50  mb 
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layer  between  25  January  and  about  3  February  and  vice  versa  thereafter 
until  the  end  of  the  analysis  period.  The  rapidity  of  the  warming  at 
25  mb  from  24  January  to  3  February  corresponded  to  the  consistent  high 
values  of  the  (ojT)  term  during  this  period.  The  occurrence,  later  on, 
of  large  scale  descent  and  rapid  temperature  increases  at  the  lower 
levels  corresponded  to  the  group  of  maximum  values  of  the  (u*T)  term 
from  about  4  February  to  9  February. 


Table  10.  Rate  of  change  of  kinetic  energy  of  the  volume  due  to  conver- 

21 

sion  from  potential  energy,  in  units  of  10  ergs  per  second. 


Date 

Time  (GOT) 

(wr) 

17  January 

0300 

2.87 

1500 

2.07 

18 

0300 

1-33 

1500 

0.93 

19 

0300 

0.05 

1500 

-1.22 

20 

0300 

-3.24 

1500 

-5.21 

21 

0300 

-5.61 

1500 

-5.  ^7 

22 

0300 

-5-93 

1500 

-6.55 

23 

0300 

-7-53 

1500 

-8.37 

24 

0300 

-8.97 

1500 

-10.34 

25 

0300 

-11.20 

1500 

-11.26 

26 

0300 

-11.68 

1500 

-11.92 

27 

0300 

-11.85 

1500 

-11.76 

28 

0300 

-11.44 

1500 

-11.26 

29 

0300 

-11.52 

1500 

-11.79 

30 

0300 

-12.17 

1500 

-12.57 
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Table  10*  Continued 


Date 

Time  (GOT) 

(uff) 

31  January 

0300 

-12*50 

1500 

-12.78 

1  February 

0300 

-12.94 

1500 

-12.59 

2 

0300 

-13.47 

1500 

-14.98 

3 

0300 

-16.09 

1500 

-18.06 

4 

0300 

-20.51 

1500 

-21.23 

5 

0300 

-21.05 

1500 

-21.87 

6 

0300 

-22.41 

1500 

-21.77 

7 

0300 

-21.49 

1500 

-21.48 

8 

0300 

-20.04 

1500 

-17.77 

9 

0300 

-15.58 

1500 

-13.31 

10 

0300 

-12.04 

1500 

-11.88 

11 

0300 

-11.54 

1500 

-IO.96 

12 

0300 

-9.96 

1500 

-9.06 

13 

0300 

-8.86 

1500 

-8.31 

14 

0300 

-7-59 

1500 

-7.57 

The  values  of  the  rate  of  dissipation  of  kinetic  energy  within 
the  volume  through  conversion  to  potential  energy,  during  the  four  sub¬ 
periods  are  shown  in  Table  11. 
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Fig.  11,  Rate  of  conversion  of  kinetic  to  potential  energy  within 
the  100-50  mb  layer  (solid  line)  and  within  the  50-25  mb 
layer  (dashed  line).  Negative  values  indicate  conversion 
from  kinetic  to  potential  energy. 
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Table  11.  Average  values  of  (afi?),  in  units  of  10  ergs  per  hour. 


Period 

(off) 

17  January  to  25  January 

-16.55 

25  January  to  4  February 

-46.53 

4  February  to  9  February 

-74.57 

9  February  to  14  February 

-37-66 

In  view  of  the  limited  volume  chosen  for  the  present  study,  no 
attempt  was  made  to  evaluate  the  total  conversion  rates  in  terms  of 
conversions  due  to  either  zonal  eddies  and  meridional  overturnings  or  in 
terms  of  conversions  in  different  wave  numbers.  The  conversion  rates 
obtained  for  the  lower  stratospheric  layer  over  the  north  American  con¬ 
tinent  should,  however,  be  viewed  in  relation  to  the  conversion  between 
available  potential  energy  and  total  kinetic  energy  over  the  rest  of  the 
Northern  Hemisphere.  In  this  connection,  the  estimates  of  conversion 
rates  of  zonal  and  eddy  available  potential  energy,  made  by  Wolfe  (i960) 
and  Reed  (1962)  are  of  particular  interest.  These  estimates  were  made 
for  the  period  25  January  to  9  February,  1957  which  is  only  a  part  of 
the  time  period  used  in  this  study.  But  the  area  chosen  by  them  covered 
the  Northern  Hemisphere  between  latitudes  30N  and  85N.  The  values  ob¬ 
tained  by  them  both  for  zonal  conversions  and  eddy  conversions  at  the 
50-mb  level  are  two  orders  of  magnitude  smaller  than  the  average  conver¬ 
sion  rates  computed  for  the  100-50  mb  layer  over  North  America,  for  the 
same  period.  In  view  of  the  differing  methods  of  averaging  used. 
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especially  in  the  computation  of  vertical  motions,  the  two  sets  of 
values  are  not  strictly  comp  sir  able.  Also,  the  area  over  which  the  in¬ 
tense  warming  occurred  was  only  a  small  portion  of  the  Northern  Hemi¬ 
sphere.  The  extremely  large  values  for  the  conversion  rates  over  the 
analysis  area  as  compared  to  those  on  a  hemispheric  scale  would,  there¬ 
fore,  indicate  that  the  rapid  conversion  to  potential  energy  in  the 
area  of  the  stratospheric  warming  was  a  major  sink  for  the  hemispheric 
kinetic  energy. 

Net  generation  of  kinetic  energy  in  the  volume 

The  values  for  the  rate  of  change  of  kinetic  energy  due  to  the 
generation  processes  (^vn)  an<l  (0o)  at  the  boundaries  of  the  volume  and 
the  conversion  process  (u?l)  within  the  volume,  are  presented  in  Table 
12,  for  immediate  reference  and  comparison  purposes. 

oh 

Table  12.  Average  values  of  the  net  generation  rate,  in  units  of  10 
ergs  per  hour. 


Period 

(0vn)+(0U)> 

(off) 

Net 

generation  rate 

17  January  to  25  January 

18.95 

-16.55 

2.40 

25  January  to  4  February 

52.97 

-46.53 

6. 44 

4  February  to  9  February 

78.19 

-74.57 

3.62 

9  February  to  14  February 

42.73 

-37.66 

5.07 

10.  ENERGY  TRANSPORT  AND  GENERATION  PROCESSES  IN  TERMS  OF  MEAN  AND 


EDDY  CONTRIBUTIONS 

In  many  previous  investigations  of  the  energy  transport  and 
transformation  processes  based  on  hemispheric  data  (e.g.,  Jensen, 
i960 )  the  usual  procedure  was  to  divide  the  terms  representing  these 
processes  into  their  mean  values  and  deviations  from  the  mean  values. 
The  latter  were  again  subdivided  to  yield  the  contributions  from  ver¬ 
tical  motions  in  a  zonal  plane,  meridional  overturnings  and  synoptic 
scale  vertical  motions  associated  with  migratory  cyclones  and  anticy¬ 
clones.  Although  a  closed  pressure  surface  around  the  earth  was  not 
used,  it  was  customary  to  neglect  the  mean  values  of  the  terms  which 
were  due  to  mean  vertical  motions  over  the  portion  of  the  hemisphere 
actually  used  in  the  analysis.  Over  an  area  much  larger  than  the 
scale  of  the  vertical  motions,  one  would  normally  expect  their  mean 
value  to  be  nearly  zero,  areas  of  downward  and  upward  motion  cancelling 
each  other.  But  the  limited  extent  of  the  analysis  area  and  depth  of 
stratospheric  layer  used  in  the  present  study,  precluded  any  such  ana¬ 
lysis  based  on  zero  mean  vertical  motion.  The  mean  u)  over  the  analysis 
area  was,  in  fact,  consistently  positive  during  most  of  the  period  and 
the  contribution  of  the  mean  motion  to  the  values  of  (ktu),  (0w)  and 
(or)  terms  could  be  quite  significant. 

Since  the  dominant  scale  of  motion  during  the  warming  was  wave 
number  two,  it  must  be  emphasized  here,  that  the  mean  motion  over  the 
analysis  area  would  only  correspond  to  what  may  be  termed  a  large  scale 
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eddy  motion  on  a  hemispheric  scale.  In  what  follows,  therefore,  the 
mean  motion  and  the  eddies  which  involve  departures  from  the  mean  motion 
refer  only  to  the  limited  portion  of  the  dominant  planetary  wave  covered 
by  the  analysis  area. 

The  procedure  adopted  for  the  treatment  of  the  (ku)),  (0u>)  and 
(oXT)  processes  in  terms  of  mean  motion  and  eddies  was  based  on  the  clas¬ 
sification  scheme  proposed  by  Benton  and  LaSeur  (1953)*  Following  the 
nomenclature  used  by  them,  the  residual  flow  that  remained  after  the  mean 
motion  was  subtracted  out  of  the  total  flow,  consisted  of  space  eddies, 
time  eddies  and  space- time  eddies.  The  space  eddies  were  defined  as 
"eddies  which  remain  after  averaging  over  the  time  period  under  consid¬ 
eration/'  Such  eddies  would,  therefore,  show  up  on  time-mean  charts. 

The  fluctuations  in  intensity  of  these  space  eddies  were  attributed  to 
other  eddies.,  Time  eddies  were  defined  as  "eddies  which  appear  in  time 
after  averaging  over  the  entire  space  field  under  consideration."  The 
time  eddies  would  reveal  any  periodicities  in  the  variations  of  the 
total  flow  pattern  during  the  particular  time  period.  Space -time 
eddies  were  defined  as  "eddies  which  disappear  when  the  motion  is  aver¬ 
aged  either  over  the  space  field  or  the  time  field  under  consideration." 
The  migratory  cyclones  and  anticyclones  of  middle  latitudes  are  examples 
of  space-time  eddies. 

For  the  purposes  of  the  present  study,  "space"  denoted  the  area 
in  Fig.  1  and  the  time  periods  used  were  the  four  subperiods  mentioned 
previously.  Accordingly,  the  expression  for  the  contributions  to  the 
vertical  transport  of  kinetic  energy  due  to  the  mean  motion  and  the 
various  eddies  of  a  two-dimensional  flow  can  be  written  as  follows. 
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(k.  .u).  .  )  =  k  .a)  +  (k  .<B  .  -  k  ,u)  +  ik.,<u,,  -  k  .to  * )  + 

^  ii  ii  J  ^  st  st  si  si  st  sty*  V  it  it  st  sty 


/ 


i  k,  ,  -  A  .  uu  .  -  A  ,  UU,  ,  <+*  A  ,  <JU  , 

y  ii  ii  si  si  it  it  st  st 

s  denotes  averaging  over  area 

t  denotes  averaging  over  time 

st  denotes  averaging  over  area  and  time 

ii  denotes  the  value  at  the  i  grid  point  at  t  =  i 

si  denotes  the  areal  mean  value  at  t  =  i 

it  denotes  the  time  mean  value  at  the  i  grid  point 

lr  (I) 

ii  ii  is  the  time  average  of  the  vertical  transport 
term,  -  S  ktu  M . ,  through  any  pressure  level. 

g  i  l 

The  first  and  second  terms  on  the  right  side  of  (37)  represent 
the  average  vertical  transport  of  kinetic  energy  through  any  pressure 
surface  due  to  the  mean  motion  and  time  eddies  respectively.  The  third 
and  fourth  terms  give  the  contributions  due  to  space  eddies  and  space- 
time  eddies.  In  view  of  the  large  scale  and  slow  changes  and  movement 
of  systems  in  the  stratosphere,  it  was  decided  to  group  these  two  terms 
together  under  the  name  "space  eddies." 

The  actual  evaluation  of  the  terms  in  (37)  for  each  pressure 
level  and  during  each  of  the  four  subperiods  was  done  as  follows. 


k  .0)  .  -  k^U), .  +  k  ,  co 


(37) 


st 


where. 


kiiu,ii 
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st  =  j  (6  hrs)  +  ^S(ka))i  AAj  (12  hrs)  + 

i-  t— 1  t=2 

•  (  5Z(ka>) ^  AA^J  (12  hrs)  +  f E(ku)) AA.  ]  (6  hrs) 

'  i t=n-l  '  • t=n 


(38) 
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where,,  N  =  number  of  hours  during  the  subperiod  and  n  is  the  number  of 
maps  during  the  subperiod.  For  example,  for  the  period  0300  GOT,  17 
January  to  0300  GCT,  25  January,  t  =  1  meant  map  for  0300  GCT  of  I7  Jan¬ 
uary,  t  =  2  meant  map  for  1500  GCT  of  17  January  and  t  =  n  meant  map  for 
0300  GCT  of  25  January. 


03 


st  na 


Euu  AA.  j  (6  hrs)  +  [ Tid .  AA  ]  (12  hrs)  + 
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•  •  +  (  5Tco .  AA 


AA^^\  (12  hrs)  +  f  Soo.  AA  J  (6  hrs)  ! 
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(39) 


Similar  expressions  could  be  written  down  for  kg>t,  0 ^  and  T 
Substitution  of  0  or  T  for  k  in  (38)  yielded  expressions  for 
(0u))gt  and  (U)T)st  for  each  pressure  level. 
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The  evaluation  of  the  space  eddy  term  given  by  the  sum  of  the 
third  and  fourth  terms  was  not  done  directly  but  obtained  as  a  residual 
out  of  the  imbalance  among  the  other  terms. 

Vertical  integration  of  expressions  (38)*  (39)*  (40)  and  similar 
expressions  involving  0  and  T,  gave  the  contributions  to  the  particular 
processes  in  the  layer  through  which  the  integration  was  carried  out. 
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The  (kto)  process  -  mean  and  eddy  contributions 

The  values  for  the  mean  term,  ^  and  the  space  and  time  eddy- 

terms  at  the  three  pressure  levels  as  well  as  their  contributions  to  the 
vertical  transport  of  kinetic  energy  through  the  layers,  are  given  in 
Tables  13  and  14  respectively. 

The  most  obvious  feature  of  these  results  is  the  relative  impor¬ 
tance  of  the  vertical  transport  of  kinetic  energy  through  mean  sinking 
motion  at  the  upper  and  lower  boundary  surfaces  of  the  layers. 

Almost  the  entire  period  of  analysis  was  characterized  by  a  net 
downward  transport  of  kinetic  energy  at  the  three  pressure  levels.  As 
indicated  in  Section  7 >  the  positive  contribution  to  kinetic  energy  in¬ 
crease  in  the  layers  as  well  as  in  the  volume,  during  the  first  two  sub¬ 
periods  wasj  in  general,  due  to  a  larger  downward  transport  rate  through 
25  mb  and  successively  lower  rates  through  the  50  and  100-mb  levels. 

The  negative  values  from  4  February  to  14  February  were  due  to  a  larger 
downward  transport  through  100  mb  and  successively  less  downward  trans¬ 
port  through  the  higher  levels.  The  values  of  a)  at  each  level,  when 
averaged  over  the  area,  were  also  positive,  except  for  two  days  in  the 
beginning  of  the  analysis  period. 

As  in  the  case  of  transport  of  kinetic  energy  through  mean  down¬ 
ward  motion,  the  transports  at  all  levels  through  space  eddies  repre¬ 
sented  by  the  covariance  of  k  and  u)  over  the  analysis  area,  were  con¬ 
sistently  downward  at  each  pressure  level.  The  positive  contribution 
from  these  eddies  for  the  period  17  January  to  25  January  was  due  to  a 
larger  downward  transport  through  the  upper  boundary  of  the  layer  than 
through  the  lower  boundary.  The  negative  contributions  arose  from 
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Table  13*  Vertical  transport  of  kinetic  energy 9  in  units  of  10  ergs 
-2  -1 

cm  sec  through  the  pressure  levels  due  to  the  mean 
motion  and  the  eddy  components.  Positive  value  indicates 
downward  transport  of  energy. 


Period 


17  January  to  25  January 


25  January  to  4  February 


4  February  to  9  February 


Level 

Mean 

motiori 

Space 

eddy 

Time 

eddy 

25  mb 

1.34 

0.77 

-0.02 

50  mb 

1.07 

0.50 

0.06 

100  mb 

0.93 

0.13 

0.03 

25  mb 

4.38 

1.12 

-0.01 

50  mb 

3.45 

1.68 

-0.08 

100  mb 

2.77 

1.90 

0.00 

25  mb 

2.40 

1.79 

0.13 

50  mb 

4.21 

3.30 

0.03 

100  mb 

6.09 

4 .66 

0.03 

25  mb 

0.28 

0.11 

0.03 

50  mb 

1.09 

0.61 

0.04 

100  mb 

2.81 

1.38 

0.01 

9  February  to  14  February 
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Table  14.  Contributions  to  the  (kw)  process,  in  units  of  erg  cm  mb 
sec  due  to  mean  motion  and  eddy  components.  Positive 
value  indicates  contribution  acted  to  increase  horizontal 
kinetic  energy. 


Period 

Layer 

Mean 

motion 

Space 

eddy 

Time 

eddy 

17  January  to  25  January 

50-25  mb 

1.07 

1.08 

-0.31 

100-50  mb 

0.30 

0.75 

0.07 

100-25  mb 

0.50 

0.84 

-0.05 

25  January  to  4  February 

50-25  mb 

3.73 

-2.22 

0.28 

100-50  mb 

1.37 

-0.44 

-0.15 

100-25  mb  -  • 

1.94 

-0.93 

-0.01 

4  February  to  9  February 

50-25  mb 

-7. 26 

-6.02 

o.4o 

100-50  mb 

-3.76 

-2.73 

-0.01 

100-25  mb 

-4.4-3 

-3.44 

0.11 

9  February  to  14  February 

50-25  mb 

-3-26 

-2.01 

-0.03 

100-50  mb 

-3-44 

-1.54 

0.06 

100-25  mb 

-3-04 

-1.53 

0.04 
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reversed  conditions  at  the  upper  and  lover  boundary  pressure  surfaces. 
The  consistently  positive  values  of  the  covariance  of  k  and  a)  at  each 
level,  also  indicated  a  preference  for  downward  motion  to  occur  in  areas 
of  high  wind  speeds,  at  least  over  the  portion  of  the  long  wave  pattern 
covered  by  the  analysis  area. 

The  relative  unimportance  of  the  time  eddies  serves  to  illus¬ 
trate  the  slow  changes  in  stratospheric  systems. 

The  (00))  process  -  mean  and  eddy  contributions 

The  values  for  the  mean  term  0^0)^  as  well  as  for  the  space  and 
time  eddies  at  each  pressure  level  are  given  in  Table  15 .  The  contribu¬ 
tions  to  the  (0m)  process  in  the  layers  are  shown  in  Table  l6« 

The  magnitudes  of  the  mean  motion  terms  show  that  large  scale 
downward  motion  over  the  analysis  area  during  the  warming  period  ac¬ 
counted  for  almost  all  of  the  generation  process  (0u)). 

The  magnitudes  of  the  contribution  from  space  eddies  were  very 
small  in  comparison  with  those  due  to  mean  sinking  motion. 

An  interesting  aspect  of  the  space  eddies  is  the  reversal  in 
sign  to  that  of  the  mean  motion  term  in  certain  layers  and  in  some  of 
the  subperiods.  Values  of  the  space  eddy  terms  at  each  pressure  level 
(Table  15)  showed  that  these  terms  were  negative  at  all  levels  that  is, 
transferred  energy  upwards,  or  increased  kinetic  energy  of  the  overlying 
layers,  from  17  January  to  4  February.  During  4  February  to  9  February, 
they  were  positive  only  at  100  mb  and  after  9  February  at  25  and  50  mb. 
The  predominantly  negative  values  of  the  space  eddy  terms  at  all  levels 
would  indicate  that  downward  motion  was  in  general  associated  with  low 
values  of  0,  that  is,  with  troughs  and  upward  motion  with  ridges. 
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reversed  conditions  at  the  upper  and  lower  boundary  pressure  surfaces. 
The  consistently  positive  values  of  the  covariance  of  k  and  u>  at  each 
level ,  also  indicated  a  preference  for  downward  motion  to  occur  in  areas 
of  high  wind  speeds,  at  least  over  the  portion  of  the  long  wave  pattern 
covered  by  the  analysis  area, 

The  relative  unimportance  of  the  time  eddies  serves  to  illus¬ 
trate  the  slow  changes  in  stratospheric  systems. 

The  (0oj)  process  -  mean  and  eddy  contributions 

The  values  for  the  mean  term  0  ^u>  ^  as  well  as  for  the  space  and 
time  eddies  at  each  pressure  level  are  given  in  Table  15 •  The  contribu¬ 
tions  to  the  (0 od)  process  in  the  layers  are  shown  in  Table  1 6. 

The  magnitudes  of  the  mean  motion  terms  show  that  large  scale 
downward  motion  over  the  analysis  area  during  the  warming  period  ac¬ 
counted  for  almost  all  of  the  generation  process  (0u>). 

The  magnitudes  of  the  contribution  from  space  eddies  were  very 
small  in  comparison  with  those  due  to  mean  sinking  motion. 

An  interesting  aspect  of  the  space  eddies  is  the  reversal  in 
sign  to  that  of  the  mean  motion  term  in  certain  layers  and  in  some  of 
the  subperiods.  Values  of  the  space  eddy  terms  at  each  pressure  level 
(Table  15 )  showed  that  these  terms  were  negative  at  all  levels  that  is, 
transferred  energy  upwards,  or  increased  kinetic  energy  of  the  overlying 
layers,  from  17  January  to  4  February.  During  4  February  to  9  February, 
they  were  positive  only  at  100  mb  and  after  9  February  at  25  and  50  mb. 
The  predominantly  negative  values  of  the  space  eddy  terms  at  all  levels 
would  indicate  that  downward  motion  was  in  general  associated  with  low 
values  of  0,  that  is,  with  troughs  and  upward  motion  with  ridges. 
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Table  15.  Vertical  transfer  of  potential  energy,  in  units  of  10  ergs 
-2  -1 

cm  sec  through  the  pressure  surfaces  due  to  the  mean 
motion  and  the  eddy  components.  Positive  value  indicates 
downward  transfer  of  energy. 


Period 


17  January  to  25  January 


25  January  to  k  February 


k  February  to  9  February 


Level 


Mean 

motion 


Space 

eddy 


Time 

eddy 


-6.22 


1.59 


25 

mb 

463.70 

50 

mb 

319.62 

100 

mb 

235.77 

25 

mb 

1262.90 

50 

mb 

1010.40 

100 

mb 

616.49 

25 

mb 

1297.20 

50 

mb 

1668.70 

100 

mb 

1691.30 

25 

mb 

435.60 

50 

mb 

84l.6o 

100 

mb 

1163.70 

-3.83 

-1.28 

-5.67 

-0.68 

-3.70 

0.50 

-5.30 

0.30 

-H.32 

0.22 

-3.90 

-1.30 

-1.40 

-0.09 

2.20 

0.03 

0.40 

0.20 

1.20 

0.00 

-3.80 

-0.06 

9  February  to  14  February 
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Table  1 6.  Contributions  to  the  (0oj)  process,  in  units  of  erg  cm  mb 
sec  due  to  mean  motion  and  eddy  components.  Positive 
value  indicates  contribution  acted  to  increase  horizontal 
kinetic  energy. 


Period 

Layer 

Mean 

motion 

Space 

eddy 

Time 

eddy 

17  January  to  25  January 

50-25  mb 

576.3 

-9.6 

11.2 

100-50  mb 

167.7 

3.7 

-1.2 

100-25  mb 

296.3 

-0.7 

2.9 

25  January  to  4  February 

50-25  mb 

1010.0 

6.4 

0.8 

100- 5Q  mb 

787.3 

12.0 

0,2 

100-25  mb 

840.3 

9.9 

0.4 

4  February  to  9  February 

50-25  mb 

-1486.0 

-10.0 

-4.8 

100-50  mb 

-45.2 

-7.2 

-0.2 

100-25  mb 

-512.3 

-7.9 

-1.7 

9  February  to  14  February 

50-25  mb 

-1624.0 

-3.2 

0.8 

100-50  mb 

-644.0 

10.0 

0.1 

100-25  mb 


-946.5 


5-5 


0.3 
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Jensen  (i960)  has  noted  the  same  type  of  relationship  in  certain  eddies 
in  the  100-50  mb  layer. 

It  was  pointed  out  in  Section  8  that  the  (0u))  process,  except 
in  the  beginning  of  the  period,  was  characterized  by  a  gain  of  kinetic 
energy  for  the  layers  through  generation  at  the  upper  boundary  surfaces 
and  a  loss  due  to  dissipation  at  the  lower  boundary  surfaces.  The  same 
argument  applied  to  the  mean  motion  contribution  to  the  (0to)  process. 

But  this  was  not  the  case,  in  general,  with  space  eddies.  The  covariance 
of  f  and  to  were  such  that,  more  often  the  layers  lost  kinetic  energy 
owing  to  dissipation  through  eddy  vertical  motions  at  the  boundary  sur¬ 
faces. 

As  in  the  case  of  the  (kio)  process,  the  contributions  due  to 
the  time  eddies  were  very  small. 

The  (off)  process  -  mean  and  eddy  contributions 

The  magnitudes  of  the  conversion  processes  within  the  layers  due 
to  mean  sinking  motions  over  the  analysis  area  and  due  to  the  eddies  are 
shown  in  Table  17. 

Compared  to  the  contributions  to  the  conversion  process  from 
the  mean  sinking  motions,  those  from  the  eddies  were  very  small.  The 
consistently  large  negative  values  of  the  conversion  process  due  to 
large  scale  sinking  motion  over  the  analysis  area  would  also  imply  the 
simultaneous  existence  of  upward  motion  over  other  relatively  cold  areas 
on  the  hemisphere. 

From  the  time  when  pronounced  warming  first  began  at  25  mb  to 
the  end  of  the  analysis  period,  that  is,  from  about  25  January  onwards, 
the  conversion  processes  due  to  the  space  eddies  were  also  acting  to 
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decrease  kinetic  energy.  In  the  few  days  prior  to  the  actual  warming, 
that  is,  during  the  first  subperiod,  the  conversions  due  to  the  space 
eddies  were  such  that  eddy  potential  energy  was  transformed  to  eddy 
kinetic  energy.  Since  only  a  portion  of  the  Northern  Hemisphere  is 
covered  by  the  analysis  area,  it  is  not  strictly  possible  to  relate  the 
conversion  processes  in  the  limited  volume  to  zonal  and  eddy  transfor¬ 
mations  on  a  hemispheric  scale. 

On  the  basis  of  the  scales  of  motion  involved,  however,  it 
would  seem  that  the  major  portion  of  the  kinetic  to  potential  energy 
conversion  occurred  on  a  scale  comparable  to  the  analysis  area.  This 
is  in  agreement  with  the  fact  that  the  dominant  scale  of  eddy  motion 
during  the  warming  was  wave  number  two  and  the  portion  of  the  hemi¬ 
sphere  under  study  was  part  of  such  a  large  scale  eddy. 
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Table  17,  Contributions  to  the  (cur)  process,  in  units  of  erg  cm’ 

mb  1  sec  1  due  to  mean  motion  and  eddy  components.  Positive 
value  indicates  contribution  acted  to  increase  horizontal 
kinetic  energy 


Period 

Layer 

Mean 

motion 

Space 

eddy 

Time 

eddy 

17  January  to  25  January 

50-25  mb 

-347.7 

4.1 

1.1 

100-50  mb 

-144.9 

1.5 

0.8 

100-25  mb 

-207.2 

2.3 

0.9 

25  January  to  4  February 

50-25  mb 

-980.6 

-4.7 

-1.1 

100-50  mb 

-431.2 

-5.9 

0.0 

100-25  mb 

-598.9 

-5.3 

-0.3 

4  February  to  9  February 

50-25  mb 

-1231+.7 

-14.6 

1.1 

100-50  mb 

-826.6 

-14.9 

0.0 

100-25  mb 

-961.9 

-14.5 

0.3 

9  February  to  14  February 

50-25  mb 

-499.9 

-1.2 

-0.6 

100-50  mb 

-496.4 

-2.6 

-0.2 

100-25  mb 

-485.2 

-2.1 

-0.3 

11.  KINETIC  ENERGY  BUDGET  OF  THE  VOLUME 


The  average  values  of  the  terms  in  equation  (12)  for  the  four 
time  periods  are  listed  in  Table  18. 

As  mentioned  previously,  the  terms  involving  transport  of  exist¬ 
ing  kinetic  energy,  namely  the  (kv^)  and  (kto)  terms,  are  negligible  in 
the  kinetic  energy  budget  of  the  volume.  This  is  also  confirmed  by  pre¬ 
vious  computations  (Jensen,  i960;  Barnes,  1962)  covering  the  Northern 
Hemisphere. 

Also,  the  time- change  of  kinetic  energy  in  the  volume  turns  out 
to  be  a  very  small  quantity  relative  to  some  of  the  other  terms.  The 
time  change  of  kinetic  energy  is  usually  neglected  in  most  investiga¬ 
tions  on  the  basis  that  it  does  not  change  significantly  over  a  period 
of  time.  In  our  case,  there  was  a  large  change  of  kinetic  energy  in 
the  volume  in  terms  of  the  initial  and  final  wind  fields.  However,  from 
the  point  of  view  of  the  physical  processes  acting  to  change  the  kinetic 
energy  in  the  volume,  the  actual  change  was  of  relatively  small  magni¬ 
tude,  Like  many  other  meteorological  quantities,  the  changes  in  kinetic 
energy  are  brought  about  by  the  small  differences  of  other  large  terms. 

During  the  first  part  of  the  period  until  4  February,  kinetic 
energy  in  the  volume  was  generated  through  vertical  motions  at  the  upper 
boundary  surface.  It  was  dissipated  through  vertical  motions  at  the 
lower  boundary  surface,  through  transfer  to  the  surrounding  atmosphere 
by  the  (0v^)  process  as  well  as  through  conversion  to  potential  energy. 

During  the  latter  part  of  the  period,  after  4  February,  the  (0u)) 
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Table  18.  Average  values  of  terms  in  the  kinetic  energy  balance  equation, 

2k 

in  units  of  10  ergs  per  hour.  Positive  value  indicates  pro¬ 
cess  acted  to  increase  kinetic  energy  of  the  volume. 


17  January- 
25  January 

25  January- 
4  February 

4  February- 
9  February 

9  February- 
14  February 

9k 

0.23 

-0.09 

-0.57 

-0.19 

(kvn) 

t 

o 

VJI 

-1.08 

0.63 

0.33 

<“">25  »b 

0.22 

0.55 

0.43 

0.04 

("“ilOO  mb 

-0.10 

-0.47 

-1.07 

-0.42 

(ktu) 

0.12 

0.08 

-0.64 

-0.38 

<**n> 

-5.91 

-14.77 

116.86 

116.55 

(^>25  mb 

147.67 

125.72 

128.85 

43.60 

(^100  mb 

-22.81 

-57-98 

-167.52 

-117.42 

(00)) 

24.86 

67.74 

-38.67 

-73.82 

(off) 

-16.55 

-46.53 

-74.57 

-37.66 

D 

-I.72 

-5-53 

-4.18 

-5.21 
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term  was  negative  while  the  (0v^)  term  became  positive  and  large.  The 
volume  gained  kinetic  energy  through  the  (0v^)  process  and  through  gen¬ 
eration  at  the  upper  boundary  surface.  It  lost  kinetic  energy  through 
dissipation  at  the  lower  boundary  surface  and  through  conversion  to 
potential  energy.  The  relatively  high  values  of  the  (0vn)  term,  also 
show  that  the  magnitudes  of  the  (00))  and  (urfT)  terms  may  not  be  represent 
ative  of  the  average  stratosphere  at  all  longitudes  during  this  period. 

The  imbalance  among  the  computed  terms  listed  in  Table  18,  pre¬ 
sumably  represents  D,  the  " frictional”  dissipation  of  kinetic  energy  In 
the  volume,  by  motions  of  smaller  scale  than  are  defined  by  the  analysis 
Estimates  of  D  are  also  shown  in  Table  18.  These  estimates  have  been 
obtained  as  differences  among  three  large  terms  and  include  to  some  ex¬ 
tent,  errors  in  the  computations.  Relative  to  the  other  processes  In 
the  volume,  the  magnitudes  of  the  dissipation  rates  appear  reasonable. 

In  view  of  the  fact  that  the  volume  covered  only  a  limited  portion  of 
the  Northern  Hemisphere,  these  rates  may  not  be  representative  of  strato 
speric  motions  on  a  hemispheric  scale. 


12,  SUMMARY  AND  CONCLUSIONS 


For  the  portion  of  the  lower  stratosphere  over  North  America, 
bounded  by  the  100  and  25-mb  pressure  surfaces,  the  terms  in  the  kinetic 
energy  balance  equation  have  been  evaluated.  Since  the  volume  under 
study  was  an  open  system,  energy  exchanges  with  the  surrounding  atmos¬ 
phere  were  taken  into  consideration.  These  exchanges  took  place  not 
only  between  the  adjacent  layers  above  and  below  the  volume  but  also 
along  the  vertical  boundaries. 

The  actual  variation  of  kinetic  energy  in  the  volume  was  found 
to  result  from  a  small  difference  among  other  large  terms.  As  far  as 
the  kinetic  energy  budget  of  the  volume  was  concerned,  terms  involving 
the  vertical  and  horizontal  transports  of  existing  kinetic  energy  were 
shown  to  be  negligible.  This  is  in  agreement  with  the  results  of  pre¬ 
vious  Investigations  on  the  advection  of  kinetic  energy  on  a  hemispheric 
scale.  Kinetic  energy  variations  in  the  lower  stratospheric  layers 
should  therefore  be  accounted  for  by  consideration  of  other  much  larger 
generation  processes.  Such  processes  consist  of  conversion  of  poten¬ 
tial  to  kinetic  energy  within  the  layers  and  generation  of  kinetic 
energy  at  the  boundaries. 

Previous  investigators  of  energy  conversions  in  stratospheric 
layers  during  other  winters,  have  differed  on  the  sign  of  the  conversion 
process,  (<jgT).  According  to  Starr  (1959)>  the  large  hydrostatic  stabil¬ 
ity  present  in  the  lower  stratosphere  indicates  that  this  region  is,  on 
the  average,  one  of  forced  vertical  motions  resulting  in  conversions  of 
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kinetic  to  potential  energy*  Jensen  (i960)  found  an  overall  positive 
transformation  from  potential  to  kinetic  energy  in  the  200-50  mb  layer, 
during  January  1958,  with  a  reverse  transformation  from  kinetic  to  po¬ 
tential  energy  in  the  100-50  mb  layer.  Barnes  (1962)  was  of  the  view 
that  the  upper  stratosphere  above  40  mb  was  characterized  by  conversion 
from  potential  to  kinetic  energy.  For  the  limited  volume  and  time  pe¬ 
riod  used  in  this  report,  the  conversion  of  energy  within  the  volume  as 
well  as  in  the  100-50  mb  and  50-25  mb  layers,  was  generally  in  the 
direction  from  kinetic  to  potential  energy. 

Further  work  with  data  for  long  enough  periods  and  at  levels 
reaching  25  mb  and  above,  is  needed  to  ascertain  the  variations  in  magni¬ 
tude  and  sign  of  the  conversion  process  in  the  stratosphere.  Previous 
investigators  of  energy  transformations  in  stratospheric  layers 
(Roberts, i960;  Hansrote  and  Lambert,  i960;  Wolfe,  i960;  Reed,  1962) 
have,  however,  agreed  on  the  importance  of  boundary  generation  processes 
in  the  energy  balance. 

'  The  (0v^)  and  (0u>)  terms  for  a  limited  volume  or  the  (00))  term 
in  the  case  of  a  hemispheric  layer  with  negligible  interactions  at  the 
equator,  represent  processes  generating  kinetic  energy  through  motions 
across  the  boundaries.  For  the  entire  atmosphere  with  zero  vertical 
motions  at  the  surface  and  at  the  top,  the  (0tu)  term  would,  of  course, 
vanish  in  the  final  summation  of  the  integral  representing  this  term. 

The  importance  of  the  (0U))  process  in  the  kinetic  energy  budget  of  tro¬ 
pospheric  layers,  has  been  demonstrated  by  Jensen  (i960).  His  evalua¬ 
tions  for  the  100- 50  mb  layer  were,  however,  based  on  the  assumption 
that  the  (0u>)  term  vanished  over  the  50-mb  surface  covering  the 
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Northern  Hemisphere. 

For  the  volume  and  period  under  study,  the  (0u))  process  at  the 
top  and  bottom  boundaries  were  such  that  the  volume  gained  kinetic 
energy  owing  to  downward  motions  at  the  25-mb  surface  and  lost  kinetic 
energy  owing  to  downward  motions  at  the  100-mb  surface.  Almost  through¬ 
out  the  period,  generation  at  the  25-mb  pressure  surface  was  an  impor¬ 
tant  and  consistent  source  of  kinetic  energy  while  dissipation  at  the 
100-mb  surface,  was  an  equally  consistent  sink  of  kinetic  energy  for 
the  volume.  The  essential  role  of  the  (0oj)  process  is  further  demon¬ 
strated  by  the  fact  that  during  the  first  part  of  the  period,  genera¬ 
tion  of  kinetic  energy  at  the  25-mb  surface,  was  the  only  important 
source  of  kinetic  energy  in  order  to  counteract  the  losses  due  to 
other  processes.  During  the  latter  part  of  the  period,  the  dissipation 
of  kinetic  energy  at  the  100-mb  surface  was  an  important  sink  in  order 
to  balance  the  large  increase  brought  about  by  generation  at  the  verti¬ 
cal  boundaries. 

The  role  of  vertical  motions  at  the  upper  boundaries  of  strato¬ 
spheric  layers,  in  the  kinetic  energy  balance,  is  not  usually  recognized 
and  is  often  neglected.  It  is  unfortunate  that  the  volume  covered  only 
a  limited  portion  of  the  Northern  Hemisphere.  It  is,  therefore,  dif¬ 
ficult  to  assess  the  importance  of  the  (0co)  process  in  relation  to  con¬ 
ditions  over  the  hemisphere.  As  pointed  out  in  Section  8,  the  large 
scale  of  downward  motions  over  North  America  was  a  significant  contribu¬ 
tory  factor  for  the  relatively  large  magnitude  of  the  (0cu)  term,  during 
the  warming  period,  Teweles  (1958)  has  described  a  simultaneous  warm¬ 
ing  that  took  place  over  eastern  Siberia  with  features  at  100  mb  very 
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similar  to  those  of  the  North  American  warming.  The  possibility  of 
another  region  over  Siberia  where  substantial  downward  transfer  of 
energy  occurred  at  the  stratospheric  levels  cannot,  therefore,  be  ex¬ 
cluded.  Both  the  areas  under  the  influence  of  the  warming  phenomena 
were  situated  north  of  latitude  50N.  The  covariance  of  0  and  o>  over 
the  hemisphere  south  of  this  latitude  could  have  been  such  that  in  the 
final  summation  over  the  hemisphere,  the  (0cd)  process  proved  to  be  neg¬ 
ligible.  Recent  computations  of  energy  conversions  over  a  major  por¬ 
tion  of  the  hemisphere  (Reed,  1962)  however,  indicate  that  such  was  not 
the  case.  The  lack  of  balance  among  the  conversion  processes  was 
attributed  to  a  "net  vertical  outflow  of  geopotential  through  the  base 
of  the  stratosphere," 

Curing  the  earlier  part  of  the  period,  horizontal  motions  at 
the  vertical  boundaries  resulted  in  a  loss  of  kinetic  energy  for  the 
volume  while  later  on,  such  motions  brought  about  an  increase  of 
kinetic  energy.  In  other  words,  the  (0vn)  term  was  a  source  of  kinetic 
energy  for  the  rest  of  the  stratosphere  during  the  period  prior  to  the 
downward  extension  of  pronounced  warming  to  the  lower  levels.  During 
the  time  when  warming  continued  throughout  the  volume  the  (^v^)  term 
acted  as  a  source  of  kinetic  energy  for  the  volume  considered  and  a 
sink  for  the  stratosphere.  The  transfer  of  energy  from  the  volume  to 
the  surrounding  stratosphere  during  the  early  stages  of  the  warming  is, 
perhaps,  not  unreasonable,  because  decrease  of  kinetic  energy  is  asso¬ 
ciated  with  the  warming  phenomenon.  What  is  surprising  is  the  large 
generation  of  kinetic  energy  during  the  period  of  pronounced  warming 
at  the  lower  levels,  due  to  horizontal  transfer  of  energy  from  the 
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surrounding  stratosphere. 

The  values  of  the  dissipation  rate  due  to  small  scale  motions , 
appear  reasonable  in  relation  to  the  magnitudes  of  the  various  genera¬ 
tion  processes.  These  dissipation  rates  were  obtained  as  differences 
among  three  large  quantities  and  include  to  some  extent,  errors  in  the 
computations.  The  relative  imporance  of  the  (0v^)  term  would  also 
indicate  that  these  dissipation  rates  might  not  be  representative  of 
the  average  stratospheric  layer  between  100  mb  and  25  mb. 

This  study  has  been  based  on  data  for  three  pressure  levels 
and  over  a  limited  area  during  a  one  month  period.  Evaluation  of  the 
results  obtained  over  North  America  in  terms  of  application  on  a  hemi¬ 
spheric  scale  is,  therefore,  a  difficult  proposition.  Moreover,  the 
time  period  chosen  for  study  is  not  necessarily  a  representative 
sample  in  time.  Nevertheless,  it  is  felt  that  the  computational  pro¬ 
cedures  adopted  for  the  present  investigation  are  fruitful.  Because 
of  the  appreciable  hydrostatic  stability  in  the  stratosphere,  the  unim¬ 
portance  of  phase  changes  of  water  and  slow  time  changes  of  temperature, 
computations  of  vertical  motion  by  the  adiabatic  method  should  yield 
values  of  at  least  the  right  sign  and  order  of  magnitude.  With  the 
increasing  availability  of  data  reaching  25  mb  and  above,  computations 
similar  to  those  used  here  can  be  extended  to  cover  larger  areas  over 
the  hemisphere.  When  using  longer  periods  of  time,  the  systematic 
effects  of  long-term  diabatic  heating  may  be  significant  enough  to  be 
incorporated  in  vertical-motion  computations. 

The  relative  importance  of  the  (0u>)  term  shows  that  in  the 
energy  balance  pertaining  to  layers  of  limited  depth,  the  values  of 


this  term  both  at  the  upper  and  lower  boundaries  should  be  considered. 
Inasmuch  as  vertical  boundaries  do  exist  for  an  atmospheric  volume  not 
covering  the  whole  earth,  the  contribution  to  the  kinetic  energy  gene¬ 
ration  due  to  the  (0v^)  term  must  also  be  estimated  in  the  absence  of 
sufficient  actual  wind  data  on  the  boundaries.  There  is  need  for 
further  studies  similar  to  the  current  one,  in  order  to  supplement  or 
modify' the  Conclusions  reached  here. 
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